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The field theoretic renormalization group (RG) and the operator product expansion are applied 
to the model of a transverse (divergence-free) vector quantity, passively advected by the "synthetic" 
turbulent flow with a finite (and not small) correlation time. The vector field is described by the 
stochastic advection-diffusion equation with the most general form of the inertial nonlinearity; it 
contains as special cases the kinematic dynamo model, linearized Navier-Stokes (NS) equation, the 
special model without the stretching term that possesses additional symmetries and has a close 
formal resemblance with the stochastic NS equation. The statistics of the advecting velocity field 
is Gaussian, with the energy spectrum E(k) oc fc 1_e and the dispersion law uj oc k~ 2+v , k being 
the momentum (wave number). The inertial- range behavior of the model is described by seven 
regimes (or universality classes) that correspond to nontrivial fixed points of the RG equations and 
exhibit anomalous scaling. The corresponding anomalous exponents are associated with the critical 
dimensions of tensor composite operators built solely of the passive vector field, which allows one to 
construct a regular perturbation expansion in e and 77; the actual calculation is performed to the first 
order (one- loop approximation), including the anisotropic sectors. Universality of the exponents, 
their (in)dependence on the forcing, effects of the large-scale anisotropy, compressibility and pressure 
are discussed. In particular, for all the scaling regimes the exponents obey a hierarchy related to 
the degree of anisotropy: the more anisotropic is the contribution of a composite operator to a 
correlation function, the faster it decays in the inertial-range. The relevance of these results for 
the real developed turbulence described by the stochastic NS equation is discussed. Key words: 
Fully developed turbulence, Anomalous scaling, Passive vector advection, Renormalization group, 
Operator product expansion. 
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I. INTRODUCTION 

It has become a commonplace to complain that theoretical understanding of turbulence remains the last unsolved 
problem of classical physics. Of course, the concept of turbulence refers to a great deal of disparate physical situations 
("almost as varied as in the realm of life," Ref. [1], p.l) and any exhaustive and ultimate "theory of turbulence," 
of course, can hardly ever be established. There is, however, a classical "list" of phenomena (or, rather, classes 
of phenomena) that represent and illustrate the main features of turbulence: existence and stability of solutions of 
hydrodynamics equations, convective turbulence, (in)stability of laminar flows and origin of turbulence, and so on. 
Those topics, which are of great practical and conceptual importance, have always remained in the focus of attention for 
theoreticians. One of them is the fully developed (homogeneous, isotropic, inertial-range) hydrodynamical turbulence. 
Detailed description of this concept and the bibliography of this old but still open subject can be found in the classical 
monographs [1-3]. 

Turbulent flows that occur in various liquids or gases at very high Reynolds numbers reveal a number of general 
aspects (cascades of energy or other conserved quantities, scaling behavior with apparently universal "anomalous 
exponents," intermittency, statistical conservation laws and so on), which support the hopes that those phenomena 
can be explained within a self-contained and internally consistent theory. Recent developments in this area are 
presented and summarized in Ref. [4]. 

The most remarkable features of developed turbulence are encoded in the single term of intermittency. This concept 
has no rigorous definition within the classical probabilistic theory; an excellent introduction can be found in Ref. [5] 
and Chap. 8 of book [1]. Roughly speaking, intermittency means that statistical properties (for example, correlation 
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or structure functions of the turbulent velocity field) are dominated by rare spatiotemporal configurations, in which 
the regions with strong turbulent activity have exotic (fractal) geometry and arc embedded into the vast regions with 
regular (laminar) flow. 

In the turbulence, such phenomenon is believed to be related to strong fluctuations of the energy flux. Therefore, it 
leads to deviations from the predictions of the celebrated Kolmogorov-Obukhov (KO) phenomenological theory [1-3]. 
Such deviations, referred to as "anomalous" or "non-dimensional" scaling, manifest themselves in singular (arguably 
power-like) dependence of correlation or structure functions on the distances and the integral (external) turbulence 
scale L. The corresponding exponents are certain nontrivial and nonlinear functions of the order of the correlation 
function, the phenomenon referred to as "multiscaling." 

Within the framework of numerous semi-heuristic models the anomalous exponents are related to statistical prop- 
erties of the local dissipation rate, the fractal (Haussdorf ) dimension of structures formed by the small-scale turbulent 
eddies, the characteristics of nontrivial structures (vortex filaments), and so on; see Refs. [1-3] for a review and further 
references. The common drawback of such models is that they are only loosely related to underlying hydrodynamical 
equations, involve arbitrary adjusting parameters and, therefore, cannot be considered to be the basis for construction 
of a systematic perturbation theory in certain small (at least formal) expansion parameter; see e.g. the remark in Ref. 
[6] . Thus serious doubts remain about the universality of anomalous exponents and the very existence of deviations 
from the KO theory. 

The term "anomalous scaling" reminds of the critical scaling in models of equilibrium phase transitions. In those, 
the field theoretic methods were successfully employed to establish the existence of self-similar (scaling) regimes and to 
construct regular perturbative calculational schemes (the famous e expansion and its relatives) for the corresponding 
exponents, scaling functions, ratios of amplitudes etc; see e.g. [7,8] and references therein. 

Here and below, by "field theoretic methods" we mean diagrammatic and functional techniques, renormalization 
theory and renormalization group, composite operators and operator algebras (operator-product or short-distance 
expansions), instanton calculus and so on. 

Of course, the analogy is far from exact. There is a big difference between the concepts of critical scaling in 
equilibrium phase transitions and anomalous scaling in turbulence. Formally speaking, in both cases one deals with 
nontrivial powers of the distance, but in the first case they are divided by the ultraviolet (UV) scale £, while in 
the second the same role is played by the integral, or infrared (IR) scale L. It was hoped that a close analogy can 
be achieved if the momentum space for turbulence be confronted with the coordinate space for critical phenomena. 
This idea was expressed in a phenomenological "dictionary," where, in particular, the viscous length I (that is, the 
UV scale of turbulence) was confronted with the correlation length (that is, the IR scale of critical phenomena), 
while the integral scale L was confronted with the molecular length; see e.g. Refs. [9]. Hence the idea of "inverse" 
renormalization group; see Refs. [10] for a recent discussion. 

The aforementioed phenomenon of multiscaling was also often opposed to critical scaling, because in the latter 
"everything is determined by just two exponents rj and vP 

It has usually been stressed that the intermittency is essentially a strongly nonlinear phenomenon, and that, 
therefore, the anomalous scaling in turbulence cannot be treated within any kind of perturbation theory. Probably 
for this reason (and because of a very low quality of some related papers) the field theoretic methods, for many years, 
have been ignored or taken with a strong skepticism by the turbulent community. The sharpened formulation of the 
state of the art, given in [1], is that the results obtained by diagrammatic methods are either wrong or can be derived 
by much simpler methods (pp. 214-215). Another, and in a sense opposite, point of view was expressed in [11]: "...the 
reason [that the problem of turbulence is still not solved] lies in the fact that the necessary field theoretic tools have 
appeared only recently." 

Although the theoretical description of the fluid turbulence on the basis of the stochastic Navier-Stokes (NS) 
equations remains essentially an open problem, considerable progress has been achieved in understanding simplified 
model systems that share some important properties with the real problem: shell models [12], stochastic Burgers 
equation [13] and passive advection by random "synthetic" velocity fields [14]. Although the shell models, discrete 
analogs of the NS equation, exhibit pronounced anomalous scaling, it has mostly been studied within numerical 
simulations. The Burgers equation with random or deterministic initial conditions has been extensively studied 
analytically, and it exhibits strong intermittency and the energy cascade. The model is interesting in itself and 
has various applications (e.g., description of the development of singularities in self-gravitating matter), but its 
relevance for the real hydrodynamical turbulence is far from obvious. In particular, Burgulcnce is a "turbulence 
without pressure" [11] and, what is more, without the energy conservation (in more than one dimensions), while the 
conservation of energy and the energy exchange between different velocity components are very important features of 
the genuine fluid turbulence. 

Probably the most important progress in the subject, achieved in the last decade of the twentieth century, was 
related to a simplified model of the fully developed turbulence, the so-called rapid-change model. The model, which 
dates back to classical studies of Batchelor, Obukhov, Kraichnan and Kazantsev, describes a scalar or vector quantity 
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(e.g. temperature, concentration of admixture particles or a weak magnetic field), passively advected by a Gaussian 
velocity field, decorrelated in time and self-similar in space (the latter property mimics some features of a real turbulent 
velocity ensemble). 

There, for the first time the existence of anomalous scaling was established on the basis of a microscopic model 
[15], and the corresponding anomalous exponents were derived within controlled approximations [16,17] and regular 
perturbation schemes [18]. Detailed review of the recent theoretical research on the passive scalar problem and more 
references can be found in [14]. 

It is important to emphasize here, that the two alternative (or complementary) analytical approaches to the rapid- 
change model are both field theoretic. In the "zero-mode approach," developed in [16,17] (see also [14]), nontrivial 
anomalous exponents are related to the zero modes (unforced solutions) of the closed exact differential equations 
satisfied by the equal-time correlation functions. From the field theoretic viewpoint, this is a realization of the well- 
known idea of self-consistent (bootstrap) equations, which involve skeleton diagrams with dressed lines and dropped 
bare terms (see e.g. Sec. 4.35 in book [8]). Owing to special features of the rapid-change models (linearity in the 
passive field and time decorrelation of the advecting field) such equations are exactly given by one- loop approximations, 
and the resulting equations in the coordinate space arc differential (and not integral or intcgro-differential as in the 
case of a general field theory). In this sense, the model is "exactly soluble." Furthermore, in contrast to the case of 
nonzero correlation time, closed equations are obtained for the equal-time correlations, which are Galilean invariant 
and, therefore, not affected by the so-called "sweeping effects" that would obscure the relevant physical interactions. 

In this connection, it should be noted that, due to the time decorrelation, in the rapid-change model there is no 
problem in relating Eulerian and Lagrangian statistics of the velocity field: they are identical. This allows one to 
perform very accurate numerical simulations in the Lagrangian frame; see [20]. 

From a more physical point of view, zero modes can be interpreted as statistical conservation laws in the dynamics 
of particle clusters [19]. The concept of statistical conservation laws appears rather general, being also confirmed by 
numerical simulations of Refs. [21,22], where the passive advection in the two-dimensional Navier-Stokes (NS) velocity 
field [21] and a shell model of a passive scalar [22] were studied. This observation is rather intriguing because in those 
models no closed equations for equal-time quantities can be derived due to the fact that the advecting velocity has a 
finite correlation time (for a passive field advected by a velocity with given statistics, closed equations can be derived 
only for different-time correlation functions, and they involve infinite diagrammatic series). 

The second systematic analytical approach to the rapid-change model, proposed in paper [18], is based on the field 
theoretic renormalization group (RG) and operator product expansion (OPE). 

To avoid possible confusion, it should be explained that in Ref. [18] and subsequent papers, the conventional 
renormalization group (and not the inverse RG in the spirit of Refs. [9,10]) was employed, which is based on the 
standard renormalization procedure (elimination of UV divergences). The solution proceeds in two main stages. In 
the first stage, the multiplicative renormalizability of the corresponding field theoretic model is demonstrated and 
the differential RG equations for its correlation functions are obtained. The asymptotic behavior of the latter on 
their UV argument {r/l) for r ^> I and any fixed (r/L) is given by IR stable fixed points of those equations. It 
involves some "scaling functions" of the IR argument {r/L), whose form is not determined by the RG equations. In 
the second stage, their behavior at r <C L is found from the OPE within the framework of the general solution of the 
RG equations. There, the crucial role is played by the critical dimensions of various composite operators, which give 
rise to an infinite family of independent scaling exponents (and hence to multiscaling) . 

Of course, these both stages (and thus the phenomenon of multiscaling) have long been known in the RG theory 
of critical behavior, where the OPE is used in the analysis of the small-(r/L) form of the scaling functions; see e.g. 
[7,8] and references therein. The distinguishing feature, specific to models of turbulence, is the existence of composite 
operators with negative critical dimensions. Such operators are termed "dangerous," because their contributions to 
the OPE diverge at {r/L) — ► 0. In the models of critical phenomena, nontrivial composite operators always have 
strictly positive dimensions, so that they only determine corrections (vanishing for {r/L) — > 0) to the leading terms 
(finite for {r/L) — > 0) in the scaling functions (the leading terms are related to the simplest operator unity with zero 
critical dimension). 

The OPE and the concept of dangerous operators in the stochastic hydrodynamics were introduced and investigated 
in detail in [23]; detailed discussion of the NS case can be found in the review paper [24], the monograph [25], and 
Chap. 6 of the book [8]. Later, the idea of negative dimensions was repeatedly introduced in connection with the 
anomalous scaling in turbulence [26], models with multifractal behavior [27] and the phenomena related to the Burgers 
equation [11,28]. 

The RG analysis of Ref. [18] has shown that dangerous operators are indeed present in the rapid-change model, and 
that their dimensions can be calculated systematically within a regular perturbation expansions, similar to the famous 
e expansion of the critical exponents. Owing to the linearity of the original stochastic equations in the passive field, 
only finite number of dangerous operators can contribute to any given structure function, which allows one to identify 
the corresponding anomalous exponent with the critical dimension of an individual composite operator. The actual 
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calculations were performed to the second [18] and third [29] orders in e (two-loop and three-loop approximations, 
respectively). Generalizations to the cases of compressible [30,31] and anisotropic [32] velocity ensembles and the 
vector advected field [33-37] have been obtained. 

The two approaches complement each other very well: the zero-mode technique allows for exact (nonperturbative) 
solutions for the anomalous exponents related to second-order correlation functions [16,38-40] (they are nontrivial 
for passive vector fields or anisotropic sectors for scalar fields), while the RG approach form the basis for systematic 
perturbative calculations of the higher-order anomalous exponents [18,29-31]. For the cases of anisotropic velocity 
ensembles or/and passively advected vector fields, where the calculations become rather involved, all the existing 
results for higher-order correlation functions were derived only by means of the RG approach and only to the leading 
order in e [32-37]. 

Besides the calculational efficiency, an important advantage of the RG approach is its relative universality: it is 
not bound to the aforementioned "solubility" of the rapid-change model and can also be applied to the case of finite 
correlation time or non-Gaussian advecting field [41-43]. 

It has been usually stressed that intermittency and anomalous scaling in turbulence are signatures of highly nonlinear 
nature of underlying dynamics. The main lesson that can be learned from the rapid-change model is, probably, that 
such phenomena can be encountered in a model, which is linear in the passive field, and in which the advecting 
velocity field is Gaussian and nonintermittent (in contrast to a more realistic case of the stochastic NS equation). 
What is more, the RG and OPE approach shows that intermittency (at least in the rapid-change model) can be 
essentially a perturbative phenomenon, in the sense that it is contained completely already in the ordinary (primitive) 
perturbation theory around a free (Gaussian) approximation. The infinite resummation of the primitive perturbation 
series, performed by the RG and OPE, gives rise to improved representations of the correlation functions, which 
reveal anomalous scaling behavior. On the other hand, these representations can be expanded back and reproduce 
the original perturbation series, no less and no more. 

Existence of exact solutions, regular perturbation schemes and accurate numerical simulations allows one to discuss, 
for the example of the rapid-change model and its relatives, the issues that are interesting within the general context of 
fully developed turbulence: universality and saturation of anomalous exponents, effects of compressibility, anisotropy 
and pressure, persistence of the large-scale anisotropy and hierarchy of anisotropic contributions, convergence prop- 
erties and nature of the e expansions, and so on. 

So far, however, aforementioned field theoretic methods have had only limited success when applied to the real fluid 
turbulence or, better to say, to the stochastic NS equation. 

The main problem of the self-consistency approach to the stochastic NS equation is the elimination of the kine- 
matic sweeping effects, which obscure relevant physical interactions and lead to a spurious strong dependence of the 
correlation functions on the integral scale. This problem was claimed solved using the so-called internal diagrammatic 
technique [44], but it leads to the violation of the translational invariance. Probably for this reason no attempts have 
been made to explicitly solve the resulting equations at least in the simplest one-loop approximation. 

The standard RG approach to the stochastic NS equation allows one to prove the independence of the inertial-rangc 
correlation functions of the viscous scale (the second Kolmogorov hypothesis) and calculate a number of representative 
constants within a regular e expansions in a reasonable agreement with experiment; sec e.g. [8,24,25] for a review and 
Refs. [45] for the most recent results. The problem of the anomalous scaling, that is, the dependence of the Galilean 
invariant correlation functions on the IR scale L, still remains open, probably due to the lack of an appropriate 
expansion parameter. Dangerous operators in that model are absent in the e expansions and can appear only at 
finite values of e. This means that they can be reliably identified only if their dimensions are derived exactly with 
the aid of Schwinger equations or Galilean symmetry. Due to the nonlinear nature of the problem, they enter the 
corresponding OPE as infinite families whose spectra of dimensions are not bounded from below, and in order to 
find their dependence on the IR scale L one has to sum up all their contributions. The needed summation of the 
most singular contributions, related to the powers of the velocity field (their critical dimensions are known exactly), 
was performed in [23] with the aid of the so-called infrared perturbation theory for the case of the different-time 
pair correlation functions. It has revealed their strong dependence on L, that physically can be explained by the 
aforementioned sweeping effects. This demonstrates that, contrary to what is sometimes claimed, these effects can be 
properly described within the RG approach, but one should combine the RG and OPE techniques and go beyond the 
plain £ expansions. Analysis of the L dependence of the Galilean invariant objects like the structure functions requires 
the explicit construction of all dangerous invariant scalar operators, exact calculation of their critical dimensions, and 
summation of their contributions in the corresponding OPE. This is clearly not a simple problem and it requires 
considerable improvement of the present techniques. 

As the intermediate step in the investigation of intermittency and anomalous scaling it is important to study 
simplified models that are, in a number of respects, closer to the real NS turbulence but still allow for analytical 
treatment. An important step is breaking the artificial assumption of the time dccorrclation of the advecting velocity 
field; sec the remarks in Refs. [21,22]. 
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In Refs. [41,42] (see also [43] for the case of compressible flow) the RG and OPE were applied to the problem of a 
passive scalar advected by a Gaussian self-similar velocity with finite (and not small) correlation time. The energy 
spectrum of the velocity in the inertial range has the form E(k) oc k x ~ 2e , while the correlation time at the momentum 
k scales as k~ 2+n . It was shown that, depending on the values of the exponents e and 77, the model reveals various 
types of inertial-range scaling regimes with nontrivial anomalous exponents, which were explicitly derived to the first 
[41,43] and second [42] orders of the double expansion in e and r\. Earlier, a similar model was proposed and studied 
in detail (using numerical simulations, in two dimensions) in [46]. Various aspects of the transport and dispersion of 
particles in random Gaussian self-similar velocity fields with finite correlation time were also studied in Refs. [47-49]. 

Another important step toward the NS turbulence is to consider the turbulent advection of passive vector fields. 
The latter can have different physical meaning: magnetic field in the Kazantsev-Kraichnan model of hydromagnetic 
turbulence in the kinematic approximation; perturbation in the linearized NS equation with prescribed statistics of 
the background field; density of an impurity with internal degrees of freedom, etc. 

Despite the obvious practical significance of these physical situations, the passive vector problem is especially 
interesting because of the insight it offers into the inertial-range behavior of the NS turbulence. Owing to the 
coupling between different components of the vector field (both by the dynamical equation and the incompressibility 
condition) and to the presence of a new stretching term in the dynamical equation, that couples the advected quantity 
to the gradient of the advecting velocity, the behavior of the passive vector field appears much richer than that of 
the scalar field: "...there is considerably more life in the large-scale transport of vector quantities," (p. 232 of Rcf. 
[1]). Indeed, passive vector fields reveal anomalous scaling already on the level of the pair correlation function [38]. 
They also develop interesting large-scale instabilities that can be interpreted as manifestation of the dynamo effect (in 
kinematic approximation); see e.g. [38,50,51]. Other important issues are mixing of composite operators, responsible 
for the anomalous scaling, and the effects of pressure on the inertial-range behavior, especially in anisotropic sectors. 

In the scalar case, the anomalous exponents for all structure functions are given by a single expression which includes 
n, the order of a function, as a parameter [16-18]. This remains true for the vector models with the stretching term 
[33,35]. In the special vector model without the stretching term, considered e.g. in [36,37], the number and the form 
of the operators entering into the relevant family depend essentially on n, and different structure functions should be 
studied separately. As a result, no general expression valid for all n exists in the model, and the anomalous exponents 
are related by finite families of composite operators rather than by individual operators [36,37]. In this respect, such 
a model is closer to the nonlinear NS equation, where the inertial-range behavior of structure functions is believed to 
be related with the Galilean-invariant operators, which form infinite families that mix heavily in renormalization; see 
[24,25]. 

Another important question that can be addressed for the passive vector model is the effects of nonlocal pressure 
terms on the anomalous scaling, in particular, the consistency of the hierarchical picture for the anisotropic anomalous 
contributions, known for the pressureless scalar [41,43,52] and magnetic [39,33] cases, with the presence of nonlocal 
terms in the closed equations for the correlation functions, caused by the pressure contributions [53] (a more detailed 
treatment is given in [37]). 

The general vector model, introduced in [35], includes as special cases the kinematic magnetic model, linearized NS 
equation and the special model without the stretching term, and thus allows one to control the pressure contribution 
and quantitatively study its effects on the anomalous scaling. The generalized model also naturally arises within the 
multiscale technique, as a result of the vertex renormalization [1]. 

Finally, it should be noted that, as the experience with the passive fields shows, the stochastic NS equation will 
show anomalous scaling already in its linearized form. Thus the results obtained for the passive vector quantity, 
with an appropriate statistics of the background field, can be considered as an approximation to the full-fledged NS 
problem, which, in principle, can be systematically improved by including nonlinear term as the perturbation. It was 
also argued that, with a proper choice of the random stirring, the passive vector model can reproduce the anomalous 
exponents of the NS velocity field exactly [54] . 

The aforementioned works, however, have mostly been confined with the Kraichnan velocity ensemble with zero 
correlation time. 

In the present paper, we consider the model of the passive vector field with the most general form of the nonlinear 
term, advected by the synthetic velocity ensemble with a finite (and not small) correlation time. The advected velocity 
field is Gaussian, with the inertial-range spectrum of the form E(k) oc k x ~ 2e and the dispersion law w oc k 2 ~ v , where k 
is the momentum (wave number). Thus we generalize the general vector model studied in [35] for the zero-correlated 
case, to the ensemble of the advecting velocity field employed, e.g., in Refs. [41-43] for the case of passive scalar. 

We shall study anomalous scaling, stability of the scaling regimes and analytically derive the anomalous exponents 
to the first order in e ~ r\. This allows us to investigate the universality of the anomalous exponents and the effects 
of compressibility, pressure, correlation time and large-scale anisotropy on the inertial-range anomalous scaling. The 
plan of the paper is as follows. 

In Sec. II we give detailed definition of the general vector model and the advecting velocity ensemble and discuss its 
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interesting special cases: the rapid-change and frozen regimes, kinematic dynamo model and linearized NS equation, 
and so on. We give the field theoretic formulation of the original stochastic problem and present the corresponding 
diagrammatic technique. In Sec. Ill we analyze the UV divergences of the model, establish its multiplicative renor- 
malizability and present the renormalization constants in the one-loop approximation. In Sec. IV we analyze possible 
scaling regimes of the model, associated with nontrivial and physically acceptable fixed points of the corresponding 
RG equations. There are seven such regimes, any one of them can be realized depending of the values of the model 
parameters (e, r\ and others). We discuss the physical meaning of these regimes (e.g., some of them correspond to 
zero, finite, or infinite correlation time of the advecting field, to the magnetic case or linearized NS equation) and their 
regions of stability in the space of the model parameters. In Sec. V we give the scaling representation for a general 
correlation function and present the general expression for the critical dimension of a composite field (operator) . Then 
we consider the most interesting case of tensor composite operators built only of the passive vector field, which play 
a crucial role in further discussion of the anomalous scaling. The critical dimension of such operator with arbitrary 
number of the fields and vector indices is presented to first order in e ~ r\. In Sec. VI we introduce the operator- 
product expansion and demonstrate its relevance to the issue of inertial-range anomalous scaling. We show that the 
anomalous exponents in our model can be identified with the critical dimensions of aforementioned composite oper- 
ators and present the leading terms of the inertial-range behavior for a number of correlation functions. In Sec. VII 
we show that the anomalous exponents of anisotropic contributions, determined by the critical dimensions of tensor 
composite fields, obey the hierarchy relations similar to those known for the passive scalar and zero-correlated cases. 
We also discuss the dependence of the anomalous exponents on the compressibility and pressure, and the effects of 
these factors on the hierarchy relations. The results obtained are briefly reviewed and discussed in the Conclusion. 



II. DESCRIPTION OF THE MODEL. THE FIELD THEORETIC FORMULATION 

Here and below, we denote x = {t, x}, d t = d/dt, di = d/dxi, and d is the (arbitrary) dimensionality of the x 
space. 

We confine ourselves to the case of transverse (divergence- free) passive vector field 8i(x), while the advecting field 
v(x) = {vi(x)} may have a longitudinal (potential) component, so that drfi = and diVi ^ 0. Thus the advection- 
diffusion equation has the form 

m + v t (1) - a v} 2) + d t r = ^0% + /*, (2.1) 

with the nonlinear terms 

V t (1) = 0, (vjBi) , vl 2) = d 3 M,) = OjdjVi. (2.2) 

Here Aq is an arbitrary parameter, V{x) is the pressure, k is the diffusivity, d 2 is the Laplace operator and fi{x) is 
a Gaussian stirring force with zero mean and correlator 

(Mx)f J (x')) = S(t - t') Ch(t/L), r = x - x'. (2.3) 

The parameter L is an integral scale related to the stirring, and CV, is a dimensionless function finite as L — > oo. 
Its precise form is unessential; for generality, it is not assumed to be isotropic. The force f%{x) maintains the steady 
state of the system and gives rise to nonzero correlation functions of the field 0. In a more realistic formulation it is 
replaced by an imposed nonzero mean value (9), see e.g. [39,33]. 

Nonlinear terms are chosen in the form of total derivatives, so that equation (2.1) is the conservation law for 9 and, 
for Aq = 1) gives the well-known equation for the magnetic field in the hydromagnetic problem. The amplitude factor 
in front of the first nonlinear term in (2.1) can be absorbed by rescaling of the velocity field and thus we set it to 

unity. The third possible structure, v/ 3 ' = di (vj6j), can be absorbed into the pressure term d{P . 

Besides the magnetic case (An = 1), the model (2.1) includes as special cases the linearized NS equation with 
prescribed statistics of the background field (Aq = — 1), and the model of passively advected vector impurity (_4 = 0), 
which possesses an additional symmetry, 9^9 + const, and has an intrinsic formal resemblance with the stochastic 
NS equation; see [36]. In these examples, the vector field has different physical interpretations: magnetic field, weak 
perturbation of the prescribed background flow, concentration or density of the impurity particles with an internal 
structure. 

Owing to the transversality condition, the pressure can be expressed as the solution of the Poisson equation, 

d 2 V={A -l)d i d j {v j e i ). (2.4) 
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For Ao = 1 (magnetic case) the pressure vanishes. In this case Eq. (2.1) also describes dynamics of the vorticity field 
advected by a given background velocity field, see e.g. [1]. 

In the real problem, the velocity v(x) satisfies the NS equation, probably with additional terms that describe the 
feedback of the advected field 9{x). We shall begin, however, with a simplified model where the statistics of v(x) is 
given: it is a Gaussian field with zero mean and correlation function 

{v i {x)v j tf))= f^J ^{p ij (k)+aQ ij (k)}^(a; ) fc)exp[-i(f-0+ik-(x-x')]. (2.5) 

Here i\j(k) = Sij — kikj/k 2 and Qij(k) = kikj/k 2 are the transverse and the longitudinal projectors, respectively, d 
is the dimensionality of the x space, a > is a free parameter (a = corresponds to the divergence-free advecting 
field, diVi = 0). For the function D v we choose 

g u 4 k*- d - 2 ^ 

D v (u), k) = — j— — f 7^12- (2.6) 

For the energy spectrum of the field v we thus obtain E(k) ~ k d ~ 1 J dcoD v (uj, k) ~ go^o fc 1_2e . Therefore, the coupling 
constant go > and the exponent e describe the equal-time velocity correlator or, equivalently, the energy spectrum, 
while the constant uo > and the exponent r/ are related to the frequency u> ~ uo K o k 2 ~ v characteristic of the mode 
k. The factor Kg in the numerator of Eq. (2.6) is explicitly isolated for the convenience later on. 

The exponents e and 77 are the analogs of the RG expansion parameter e = 4—d in the theory of critical behavior, and 
we shall use the traditional term "e expansion" for the double expansion in the e-r/ plane around the origin s = r\ = 0, 
with the additional convention that i] — 0(e). The IR regularization is provided by the cut-off in the integral (2.5) 
from below at k ~ to, where m ~ 1/L is the reciprocal of the integral scale. Dimensionality considerations show that 
the coupling constants go, u are related to the characteristic UV momentum scale A ~ l/£ by 

. 90 ~A 2£ , w -A". (2.7) 

The model (2.5), (2.6) contains two special cases that possess some interest on their own. In the limit u — > 00, 
g'o = go /uo = const we arrive at the rapid-change model: 

D v (LO,k)^g' K k- d -^ ( = 2e-rj, (2.8) 

and the limit uo — > 0, go = const corresponds to the case of a "frozen" velocity field (or "quenched disorder"): 

D v {aj,k)^g i4k- d+2 - 2s irS{uj), (2.9) 

then the velocity correlator (2.5) is independent of the time variable t — t' in the t representation. 

The stochastic problem (2.1), (2.3), (2.5) is equivalent to the field theoretic model of the extended set of three fields 
$ = {6',6,-v} with action functional 

S*($) = e'D f 9'/2 + 9' \-d t + Kod 2 9 - V {1) + A V {2 A - vL»- 1 v/2 (2.10) 

with V {1 ^ from (2.2). This means that statistical averages of random quantities in the original stochastic problem 
can be represented as functional averages with the weight exp5($). The first five terms in Eq. (2.10) represent the 
so-called Martin-Siggia-Rose action for the stochastic problem (2.1), (2.3) at fixed v (see, e.g., [8,24,25] and references 
therein), while the last term represents the Gaussian averaging over v. Here Df and D v are the correlation functions 
(2.3) and (2.5), respectively, 9' = 9' i (x) is an auxiliary transverse vector field, the required integrations over x = (t,x) 
and summations over the vector indices are implied, for example, 

9'd t 9= f dtdx9i(t,x)d t 9i(t,x). 

The pressure term can be omitted in the functional (2.10) owing to the transversality of the auxiliary field: 

J dx9%V = - J dxVd l 9' l = 0. 

Of course, this does not mean that the pressure contribution can simply be neglected: the field 9' acts as the transverse 
projector and selects the transverse part of the expression in the square brackets in Eq. (2.10). 
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The model (2.10) corresponds to a standard Feynman diagrammatic technique with the triple vertex 9'[— + 
AoV^ 2 '] and bare propagators in the frequency-momentum (w-k) representation 

(^>0= (_i5+^2 ) . M-)o= (^+^4)' W>0= ' (2- 11 ) 

where Cy(k) is the Fourier transform of the function Cy (r/L) from (2.3); the bare propagator (u,Uj)o is given by 
Eq. (2.5). 



III. UV RENORMALIZATION. RG FUNCTIONS AND RG EQUATIONS 

The analysis of UV divergences is based on the analysis of canonical dimensions. Dynamical models of the type 
(2.10), in contrast to static models, have two scales, i.e., the canonical dimension of some quantity F (a field or a 
parameter in the action functional) is described by two numbers, the momentum dimension d F and the frequency 
dimension d F ; see e.g. [8,24,25]. They are determined so that [F] <~ [_L] _d F[T] _d F, where L is the length scale and T 
is the time scale. The dimensions are found from the obvious normalization conditions d\. — — c?£ = 1, d% = d% = 0, 
dt = d-t — 0, d% = — df = 1, and from the requirement that each term of the action functional be dimensionless (with 
respect to the momentum and frequency dimensions separately). Then, based on dp and d F , one can introduce the 
total canonical dimension dp — d F + 2d F (in the free theory, 8t ex 8 2 ), which plays in the theory of rcnormalization 
of dynamical models the same role as the conventional (momentum) dimension does in static problems. 

The dimensions for the model (2.10) are given in Table I, including the parameters which will be introduced later 
on. From the Table it follows that the model is logarithmic (the coupling constants go, uq are dimensionless) at 
e = 1] = 0, so that the UV divergences in the correlation functions have the form of the poles in e, 77 and their linear 
combinations. 

The total canonical dimension of an arbitrary 1-irreducible correlation function T ~ ($ • • • < &)i-i r is given by the 
relation d r = d T + 2d f = d + 2 - AM*, where 7V$ = {N 9l N e >,N v } are the numbers of corresponding fields entering 
into the function T, and the summation over all types of the fields is implied. The total dimension dr is the formal 
index of the UV divergence. Superficial UV divergences, whose removal requires counterterms, can be present only 
in those functions T for which rfr is a non-negative integer. 

Analysis of the divergences in our model can be augmented by the following considerations: 

(i) From the explicit form of the vertex and bare propagators it follows that Ng> — Ng = 2N for any 1-irreduciblc 
correlation function, where iVo > is the total number of bare propagators {99) entering into the function. Therefore, 
the difference Ng> — Ng is an even non-negative integer for any nonvanishing function; cf. Refs. [18,41,43]. 

(ii) For any model with the Martin-Siggia-Rose-type action, all the 1-irreducible functions with Ng> = contain 
closed contours of retarded propagators {99') and vanish; see e.g. [8,24,25]. 

(iii) If for some reason a number of external momenta occurs as an overall factor in all the diagrams of a given Green 
function, the real index of divergence d' T is smaller than dr by the corresponding number (the correlation function 
requires counterterms only if d' r is a non-negative integer); see e.g. [8,24,25]. In the model (2.10), the derivative 8 at 
the vertex 9'[-V^ + A V^} can be moved onto the field 9' using the integration by parts, which decreases the real 
index of divergence: d' r = rfr — Ng/, and the field 9' enters the counterterms only in the form of the derivative 89'. 

From the dimensions in Table I we find d r = d + 2 - 7V V + Ng - (d + l)N e , and d' T = [d + 2)(1 - N g >) + Ng - N v . 
Bearing in mind that Ng> > Ng we conclude that for any d, superficial divergences can exist only in the 1-irreducible 
functions (0')i_ ir with d r = d' r = 1, (0'0)i_ ir with d r = 2, d' r = 1 and (0'0v)i_ ir with d r = 1, d' T = 0. The 
corresponding counterterms necessarily reduce to the forms 89' (which vanishes identically), 9'8 2 9, 9'V^ and 9'V^ 
with U^ 1 ' 2 ) from (2.2). The structure 9'8t9 does not contain a spatial derivative, while 9'V^ with = 8i {vj9j) 
has the form of a total derivative and vanishes after the integration over x. 

We thus conclude that our model (2.10) is multiplicatively renormalizable and the corresponding rcnormalized 
action has the form 

= 9'D f 9'/2 + 9' \-d t + kZ k 8 2 9 - Z X V { ^ + Z 2 AV (2) ] - v£>- 1 v/2. (3.1) 

Here and below g, u, k and A (without a subscript) denote the renormalizcd analogs of the corresponding bare pa- 
rameters (with the subscript 0). The correlation function D v in (3.1) should be expressed in terms of the renormalized 
parameters and Zi = Zi(g, u, A, a, e, rj, d) are the renormalization constants. The introduction of the counterterms is 
reproduced by the multiplicative renormalization of the velocity field, v — > Ziv, and the parameters g , it , n a and 
Ao in the action functional (2.10): 
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kq = kZ k , u a = u/jP Z u g = gn E Z g , A = Z A A. (3.2) 

Here /j, is the reference mass (additional arbitrary parameter of the renormalizcd theory) and the renormalization 
constants in Eqs. (3.1) and (3.2) are related as follows: 

Z g = ZjZ~ 2 , Z U = Z~ 1 , Za = Z 2 Z 1 1 . (3.3) 

The first two relations in Eq. (3.3) result from the absence of the renormalization of the term with D v in (3.1). No 
renormalization of the fields 6, 9' and the parameters m ~ 1/L and a is required, i.e., Zg = 1 and so on. 

We have calculated all the renormalization constants in the one-loop approximation (first order in g). The resulting 
expressions are rather cumbersome: they are given by infinite series in the parameter u with the terms containing the 
poles in 2e + sr] with s = 1, 2, . . ., cf. Refs. [41,43] for the scalar case. For this reason, below we give them only for the 
special case r\ = and arbitrary e. It is important here that the parameter e alone provides the UV regularization 
for the theory, so that the constants Z remain finite at 77 = 0. In the minimal subtraction (MS) scheme they have the 
form: 



where 



Zk = 1 ~ 4(l + u)e {X + {a 1)y} + ° (ff2) ' (3 ' 4C) 

X 1 ( 1 -- 4 ) 2 2 1 2( U + 2)(1-^ 

d d(u+l) + (« + l)d(d + 2) ' [ ' 

= l + A-A 2 (I -A) 2 2 2(1 -A) 

y d d(d + 2) d(u + l) d(d + 2)(u+l)' K ' ' 

Here and below S d = S d /(2n) d and S d = 2-K d ' 2 /T(d/2) is the surface area of the unit sphere in d-dimensional space. 

The explicit expressions (3.4) illustrate some general properties of the renormalization constants ^1,2, valid to all 
orders in g. 

For the magnetic case A — 1 we have Z\ — Z 2 and, therefore, Zj± — 1 and the parameter Aq is not renormalizcd: 
„4o = A = \. This is a consequence of the relation di — = for the vectors (2.2), which expresses the 

transversality of the vertex 9[ v} 1 ^ ~ with respect to the index of the field 9'. The same property holds for 

all the diagrams of the 1-irreducible function {9'9v)\~- u and, as a result, for the corresponding countertcrm, cf. Refs. 
[55,56] for the case of the active magnetic field interacting with the NS velocity field. 

In the rapid-change limit (u — > 00, g/u 2 = const) we obtain Z\ = Z 2 = 1 due to the fact that all the diagrams of the 
function (#'0v)i_i r contain effectively closed circuits of retarded propagators {99') and therefore vanish; it is crucial 
here that the correlation function (2.8) is proportional to the 5 function in time representation. On the contrary, in 
the frozen limit (it — > 0, g/u — const) the constants Zi_ 2 remain nontrivial. 

One also obtains Z\ = Z 2 = 1 for A = a = 0. In this case, the derivative d at the only vertex 9'jV^ = 9' i (vd)9 i 
can be moved onto either fields 9 and 8', so that the real index of divergence takes on the form dp = dp — Ng — Ngi 
(we recall that d' r = dp — Ng> for general A and a). This gives d' r = — 1 for Ng = Ng> = N v = 1, so that the function 
(8'6v)i-i T is UV finite; cf. Ref. [41] for the scalar case. In other words, the counterterm to the vertex must include 
two derivatives (one for the field 9 and one for 9'), which is forbidden by the dimensionality considerations. 

Finally, from Eqs. (3.4) it follows that Z\_ 2 = 1 for A = 1 and a = 0. We found no general explanation for 
this relation, but checked that it remains true in the two-loop approximation, so that we can only guarantee that 
Z h2 = l + 0(g 3 ). 

Let W(eo) be some correlation function in the original model (2.10) and Wn(e,ii) its analog in the renormalized 
theory with action (3.1). Here e is the complete set of bare parameters, and e is the set of their renormalized 
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counterparts. In the following, we shall not be interested in the correlation functions involving the velocity field 
v. Then the relation S(8, 8', Z 2 v, e ) = Sr(6,8' , v, e, /tt) for the action functional yields W(e ) — Wn(e,n) for any 
correlation function of the fields 6', 9; the only difference is in the choice of variables and in the form of perturbation 
theory (in g instead of go). 

We use to denote the differential operation /id^ for fixed eo and operate on both sides of this equation with it. 
This gives the basic RG equation: 

V RG W R {e,n)=Q, (3.6) 

where T>rq is the operation T>^ expressed in the rcnormalized variables: 

D BG =V„ + [3 g d g + f3 u du + PaOa - lJ> K . (3.7) 

In Eq. (3.7), we have written V x = xd x for any variable x, the RG functions (the (3 functions and the anomalous 
dimensions 7) arc defined as 

7i = X> M lnZi (3.8) 

for any renormalization constant Zi and 

p g = V ft g = 2g[-e + j K -j 1 ], (3 U = V^u = u[-r] + j K ], p A = V^A = A[ 11 ~ l2 ]- (3.9) 

The relations between (3 and 7 in (3.9) result from the definitions and the relations (3.3). 

For the basis anomalous dimensions from the definitions and the expressions (3.4) in the one-loop approximation 
we obtain: 

-gS d f A(l-A) (1-^)1 „, 2n /o 1 n^ 

* = 2dirw r + W - a W } + 0{9 } ' (3 - 10a) 



7, - U - fe4 + \ + 0(g% (3.10b) 



2d(l + u) 2 \ {d + 2) A(d + 2) 
gS, 



1* 



2(1 + u) 



[X + (a - l)y} + 0(g 2 ) (3.10c) 



with X and y from (3.5). It is also worth noting that the knowledge of the constants Z at r\ — is in fact sufficient to 
calculate the functions (3.9) for all rj, e because the anomalous dimensions (3.10) are independent of 77, e, at least 
in the one-loop approximation. This fact essentially simplified the two-loop calculation of the anomalous dimensions 
for the scalar case performed in Ref. [42] . 



IV. FIXED POINTS AND SCALING REGIMES 



It is well known that possible scaling regimes of a renormalizable model are associated with IR attractive fixed points 
of the corresponding RG equation; see e.g. Ref. [8] . Roughly speaking, in solving the RG equation all the renormalized 
coupling constants gt (that is, all dimensionless parameters of the model) are replaced by the corresponding invariant 
charges g~i(s), where s = k/fj, in the momentum representation or s = 1/fJ.r in the coordinate representation. The 
invariant charges are determined as the solutions of the following Cauchy problem 

V s gi(s) = &({&(*)}), 9i0-)=9i for all*. (4.1) 

Here §i(s) is the full set of the invariant charges and /3j = V^gi are the corresponding (3 functions. In the IR asymptotic 
range (s — ► 0) the invariant charges tend to IR attractive fixed points of the system (4.1). The coordinates g« of 
the fixed points are found from the equations (3i{{gj*}) = for all i. The type of a fixed point is determined by the 
matrix £1 = {fijj = dfli/dgj}: for an IR attractive fixed point the matrix is positive definite, i.e., the real parts of 
all its eigenvalues are positive. 

In our case, the coordinates 5*, u* , A* of the fixed points are found from the equations 
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fig (.9* 1 ^* 5 ) fill ($* 5 ^* ) *^-* ) ft A. G?* 5 ^* •) ) 



(4.2) 



with the (3 functions given in Eq. (3.9). The coordinates of the fixed points and the elements of the corresponding 
matrices f2 depend on the remaining free parameters: e, rj, d and a. 1 

Below we list all possible fixed points of the system (4.2), giving their coordinates in the one-loop approximation, 
that is, to first order in e and r\. We shall also present the inequalities that determine the regions (in the space of 
parameters e, 77, d and a) where those points are IR attractive. It should also be kept in mind that admissible fixed 
points should satisfy the relations g* > 0, > 0, which follow from the physical meaning ot these parameters (g is 
the amplitude of a pair correlation function and u is the ratio of the diffusivity and viscosity coefficients). 

First of all, the trivial fixed point 

.9* = ?i* = 0, A* arbitrary (4.3a) 

should be mentioned. The corresponding matrix O is diagonal with the diagonal elements (eigenvalues) 

Ai =0, A 2 = -2e, A 3 = -77, (4.3b) 

so that the point (4.3a) is IR attractive for e < 0, 77 < 0. Since for g = u = all the three (3 functions (3.9) vanish 
simultaneously, the value of A at this fixed point remains arbitrary. This degeneracy is reflected in the vanishing of 
Ai. 

We shall discuss the physical meaning of the point (4.3a) later and now turn to nontrivial fixed points. Their 
analysis is simplified by the observation that the one-loop function Pa factorizes into a part that depends only on A 
and a part that depends only on g, u: 

PA = 1 £^ { A*-lXA-a) + Otf), (4.4) 

see Eq. (3.9), (3.10). Thus all possible values of A* are found from the equation (3a — regardless of the values of 
the other couplings, provided g is nonzero and g, u are not infinite: 

A, = -1, 1, a. (4.5) 

Furthermore, the elements d(3A/dg\A=A, = d&Ai i ® u \a=A, = of the matrix vanish for these values of A* regardless 
of the values of g, u. Thus the matrix Q is block-triangular and one of its eigenvalues coincides with the diagonal 
element ft a = 9/3/i/cVl. Its sign is completely determined by the value of A*, so that we can check a necessary 
condition VIa > for a fixed point to be IR attractive regardless of the values of g, u: it is satisfied for A* = —1 for 
all a, for A* = 1 if a < 1, and for A* = a if a > 1. This is easy to understand geometrically: ft a > for the leftmost 
and rightmost points and FIa < for the point lying in between them. 

Therefore, for any a the condition flA > is simultaneously satisfied by two fixed points: A* = —1 and A* = 
max {1, a}. Factorization of the function (4.4) allows one to analyze the solution of the equation (4.1) for the 
invariant charge A(s) independently of the remaining equations for the other invariant charges: A(s) will be attracted 
by the leftmost fixed point A* = —1 if and only if the initial condition A = A(l) lies to the left of the unstable fixed 
point, A < min {1, en}, and by the rightmost point A* — max {1, a} if and only if A > min {1, a}. 

Besides the finite values (4.5), there is one more possibility that formally corresponds to A* = 00. More accurately 
it can be revealed by the change of variables a = 1/A, y = gA 2 ; then a* = and is finite. It describes the situation 
when the action (2.10) contains the only vertex 9'V^ 2 \ Since in this vertex the derivative can be moved onto either of 
the fields 9' and v [see Eq. (2.2)], such a model is multiplicatively renormalizable (the vertex O'V^ is not generated 
by the renormalization) and Z2 = 1 identically. Thus the result a* = is exact to all orders. However, one can easily 
check that the eigenvalue of the matrix Q, equal to the diagonal element Q a = df3 a /da, is always negative, so that 
this point cannot be IR attractive. We shall not discuss it in what follows. 

To conclude the analysis of the equation (3a — 0, it remains to note that the result A* = 1 in Eq. (4.5) is exact (a 
consequence of the relation Z\ = for A = 1, see Sec. Ill), while the other two can have corrections of order e ~ r/ 



1 Formally, a can be treated as the fourth coupling constant. The corresponding /3 function /3 a = V^a vanishes identically 
owing to the fact that a is not renormalized. Therefore, the equation j3 a = gives no additional constraint on the values of 
the coupling constants at a fixed point. 
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and higher. The only exception is the result A* = for a = 0, which is also exact due to the relation Z\ = 1 for 
a = 0. 

Substituting the values (4.5) into the functions f3 g , (3 U from (3.9) and solving the equations f3 g — f3 u — gives the 
values of the remaining coordinates g* and w*. For each value of .4*, there are two solutions. For the first of them, 
it* = 0. This result is exact to all orders of the e expansion since the function (3 U for u = vanishes identically, see 
Eq. (3.9). Substituting the value = to f3 g and solving the equation (i g = gives g*. We shall denote such fixed 
points by Q, as they correspond to the case of "quenched disorder," see the comments to Eq. (2.9). For the second 
variant, g* and are both nonzero; we shall denote such fixed points by F ("finite" correlation time of the velocity 
field). 

Thus we arrive at six nontrivial fixed points Q~ , Q + , Q a , F~ , F + , F a , where the superscripts correspond to the 
values of A* in Eq. (4.5). Below we give the coordinates of these points in the one-loop approximation, the eigenvalues 
Ai.2,3 of the matrix £1 (Ai always corresponds to the diagonal clement df3_A/dA), and the inequalities that determine 
the regions where the points are admissible (IR attractive and satisfy > 0, u* > 0): 



with the eigenvalues 



2de 

Q + : = - _ - , u* = 0, A* = 1 (4.6a) 



2e(l - a) 



Al (d-l)(d + 2)' 

A 2 = 2e, 

A 3 = (d-l-«)g~^-l) > (4 . 6b) 

admissible for 

a < 1, (d - 1 - a)e > {d - l)r]. (4.6c) 



2d(d + 2)e n „ , 

Q ~-- 9 *= S d{ d-2a)(d-l) > ^ = °' = ^ 



with the eigenvalues 

2e(a + l) 



Al_ (d-l)(d-2a)' 
A 2 = 2e, 

[-2 + d-d 2 +a(3d-2)]e-(2a-d)(d-l)?7 
A3 = (2a-d)(d-l) ' (47b) 

admissible for 

e>0, a<d/2, (d - 2a) (d - 1)77 < [d 2 - d + 2 - a(3d - 2)]e. (4.7c) 



Q Q : .9, = — ^tw^, 2V = 0, A = a (4.8a) 

d^ - 3 + a(d + 1)(1 + a — a 2 ) 



with the eigenvalues 

(a 2 - l)e 



1 ~ d 2 + a(d + 1)(1 + a - a 2 ) - 3 ' 
A 2 = 2e, 

_ [-3 - a + d 2 + a 2 (l + d)(l - a)]g - [-3 + d 2 + a(l + d)(l + a - a 2 )]// 
3_ d 2 +a(l + d)(l + a-a 2 )-3 ' [ ' 



admissible for 
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e > 0, a > 1, d 2 - 3 + a(d+ 1)(1 + a - a 2 ) > 0, 
d 2 - 3 - a + a 2 (d+ 1)(1 - a) 



d? -3 + a(d+l)(l + a- a 2 ) 



with the eigenvalues 



where 



2ad(?? - 2e) 2 
(a + d- l) 2 (e -ry)' 



(l + a-d)e+{d-l)r) 
(d + a-l)(e-ri) ' 



Ai 



A 2 = 



2(1 -cQ(g- 77) 
a(d+2) ; 



Wi - y/W? - Wa 
2a(2e - 77) 



Wi + y/w? - W 2 

2a(2e - r]) 



Wi = (2 + 6a - 2d)e 2 + 5(d - 1 - a)er/ - 3(d - l)ri 2 , 

VK 2 = 8a(2e - r?)((l + a - d)e + (d - l)r/)(2e 2 - 3e?7 + rj 2 ), 



admissible for 



d — a — 1 

£>0, r/ > 0, a < 1, e > 7] > — - — - — e, 
(2 + 6a - 2d)e 2 + 5(d - 1 - a)e?7 - 3(d - l)ry 2 > 0. 



2d(2 + d)(ad-2)(ri-2e) 2 
(l-a + d) 2 (d-2) 2 (e-r?)' 



with the eigenvalues 



where 



admissible for 



(2 + a{2 - 3d) - d + d 2 )e + (2a -d)(d- l)rj 
(a-d-l)(d-2)(e-r/) ' 

2(1 + a)(e-rj) 



Ai 

A 2 = 
A 3 = 



-2 + ad 



~ \/W? - W2 
2(-2 + ad)(2e-r;)' 

y^i + yw? - w 2 

2(-2 + ad)(2e - r7) ' 



W t = 2(-6 + d-d 2 + a(-2 + 5d))e 2 

+ 5(2 + a(2 - 3d) - d + d 2 )er/ + 3(2a - d)(-l + d)r/ 2 , 
VK 2 = 8(-2 + ad)(2e - ri) 2 (e - 77) 

x ((-2 + d - d 2 + a(-2 + 3d))e - (2a - d)(-l + d)r}), 



£>r], 2/d<a<(d+l), 2e>r), W 1 > 0, 
(2 + a(2 - 3d) + d(d - l))e < (d - l)(d - 2a)??. 
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2ad(2 + d) 2 (r)-2e) 2 

9 * ~ {d 2 - 3 + a 2 (l + d)(l - a) + a(2d + 3)) 2 (e - r,) ' * ~ 
(3 + a - d 2 - a 2 (l + d)(l - a))e + (d 2 - 3 + a(l + d){l + a - a 2 ))r] 



with the eigenvalues 



{d 2 - 3 + a 2 (l + d)(l -a)+ a{2d + 3))(e - rj) 
Xi= (* 2 -l)(e-v) 



(4.11a) 



where 



a(2 + d) ' 

= w x - ^W 2 - w 2 

2 2a(2 + rf)(2e-?7)' 
Wi + yJW 2 - W 2 

- 2a(2 + d)( 2e -^) ' (4 - llb) 



W 1 = 2(3 - d 2 - a 2 (l + d) + a 3 (l + d) + a(5 + 2d))e 2 

- 5(3 + a - d 2 - a 2 (l + d) + a 3 (l + d))sr] + 3(3 - d 2 - a(l + + a - a 2 ))?? 2 , 
VK 2 = 8a(2 + rf)(2e-?7) 2 (e-77) 

x ((3 + a - d 2 - a 2 {l + d)(l - a))e 

+ (-3 + d 2 + a(l + + a - a 2 )r/), (4.11c) 



admissible for 



a > 1, £ > jj, 2e > ?7, Wi > 0, 

(d 2 - 3 + a 2 (l + d)(l - a) + a(2d + 3)) > 0, 

(3 + a - d 2 - a 2 (l + d)(l - a))e + (d 2 - 3 + a(l + + a - a 2 ))?] > 0. (4.11d) 

However, the above list is not exhaustive: besides the fixed points with non-infinite values of <?* and u*, there are 
points for which these parameters tend to infinity. They can be revealed by the change of variables x = g/u, w = 1/u. 
The corresponding (3 functions are obtained by the chain rule: 

X = V^x - {l/u)P g - (g/u 2 )p u , t3 w =V^w = -{l/u 2 )f3 u , (4.12) 

with (ig^ u from (3.9) and the anomalous dimensions (3.10) expressed in the variables x, w. Solving the equations 
fix = Pw = Pa = gives three fixed points with finite values of £*, w*, which simply express the points F ± ' a in the 
new variables. Besides them, there are two fixed points with w* = 0, left out in the analysis performed in terms of g, 
u. It is clear from (2.8) that the choice w = corresponds to the rapid-change limit of our model. The dimension j K 
from Eq. (3.10) remains finite for w — 0, while 71 and 72 vanish. In fact, they vanish to all orders of the perturbation 
theory in x oc g owing to the exact relation Z\ = Z 2 = 1 that holds in the limit (2.8); see the discussion in Sec. III. As 
a result, the function (3a vanishes identically for w = and the coordinate A* remains arbitrary at such fixed points; 
see the remark below Eq. (4.3a). 



There is a trivial fixed point 



with the eigenvalues 



x* = w* = 0, A* arbitrary (4.13a) 



Ai = 0, A 2 = 77, A 3 = n - 2s, (4.13b) 
and a non-trivial fixed point, which we shall denote by R A in what follows, 

2d(d + 2)(2s-ri) 



R . x # 



A(l + a)d + d 2 + a(3 + d) - A 2 (-l + a + ad) - 3 ' 
w* = 0, A* arbitrary, (4.14a) 



with the eigenvalues 
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Ai=0, A 2 = -2( £ -r ? ), A 3 = 2e-77, 



(4.14b) 



admissible for 

r]>e>r]/2, A(l + a)d + d 2 + a(3 + d) - A 2 (-l + a + ad) - 3 > 0. (4.14c) 

The notation R A implies that this point corresponds to the "rapid-change" regime and that A* = A remains a free 
parameter. 

Triviality of the points (4.3a) and (4.13a) implies the absence of anomalous scaling; they correspond to diffusive- type 
regimes, for which the convection (that is, the nonlinearity in Eq. (2.1)) can be treated within ordinary perturbation 
theory and the standard methods of the homogenization theory apply. More detailed discussion of such fixed points 
(in particular, the difference between the "quenched" and "rapid-change" trivial fixed points) can be found Ref. [41] 
for the example of the passive scalar field, advected by the velocity ensemble (2.5). 

On the contrary, the nontrivial fixed points Q^' a , F ±,a and R A describe non-diffusive asymptotic regimes, in 
which the competition of the diffusive and convective terms in Eq. (2.1) produces anomalous scaling behavior. The 
corresponding anomalous exponents will be presented in the next Section, and now we shall discuss the interplay 
between the possible scaling regimes. 

Seven nontrivial IR attractive fixed points correspond to seven possible scaling regimes; only one of them can be 
realized when the values of all parameters a, e, r/, d and A are given, regardless of the values of the amplitudes 
g and u (see below). In this sense, the asymptotic behavior in our model is universal, and (as we shall see) the 
anomalous exponents depend only on the values of the exponents e and r\ in the velocity correlation function and on 
the parameters a (for all regimes) and A (for R A ), but they do not depend on the coupling constants g and u. 

Indeed, let us fix the value of A. Then the solution of the RG equation (3.6) can be attracted by either of the 
following three sets: by the set {F~, Q~ , R A } if A < min {1, a}, by the set {F + , Q + , R A }, if a < 1 and A > a, and 
by the set {F a ,Q a ,R A } if a > 1 and A > 1. (In all these cases, the value of A* for the fixed point R A in Eq. (4.14a) 
simply equals to A). 

For any of these three situations, only one fixed point in the list {F, Q,R} can be IR attractive for given e and 
rj. Indeed, the analysis of the inequalities that determine the admissibility regions for these points shows that these 
regions adjoin each other without overlaps or gaps. The common boundary of the admissibility regions for the points 
F and R is e = rj [one of the admissibility conditions for R A is e < i], see Eq. (4.14c), while one of the admissibility 
conditions for any of the points F ±a is e > rj, see Eqs. (4.9d), (4.10d) and (4.11d)]. The common boundary of the 
admissibility regions for the points F and Q is (d — 1 — a)e = (d — for the pair Q + , F + [see Eqs. (4.6c) and 
(4.9d)], (d 2 - 3 + a(d + a~ a 2 ))r] = (d 2 - 3 - a + a 2 (d + 1)(1 - a))e for the pair Q~ , F~ [see Eqs. (4.7c) 
and (4.10d)] and (d 2 - 3 + a 2 (l + d)(l - a) + a(2d + 3)) = for the pair F a , Q a [see Eqs. (4.8c) and (4.11d)]. 

Therefore, for any given set of parameters a, e, 77, d and A, the RG trajectory can be attracted by the only one 
possible nontrivial fixed point, regardless of the values of the parameters g, u [that is, the initial data for the Cauchy 
problem (4.1)]. If the trajectory is attracted by a trivial point, the model will show diffusive-type behavior. Finally, 
if there is no IR attractive fixed point for the given set of parameters, no definitive conclusion can be made about the 
asymptotic behavior of the model within the framework of the e expansion. In particular, if the trajectory comes to 
a region where the parameters g, u are negative, one can think that the steady state of our system becomes unstable 
(by analogy with the RG theory of critical phenomena, where such behavior is usually interpreted as a first-order 
phase transition). 

The admissibility regions in the e--q plane are shown in Figs. 1 for a number of values of the parameters a and d. 
In the one-loop approximation, their boundaries are always given by straight rays starting at the origin s-r/. Some of 
them, however, can be affected by the higher-order corrections, so that the gaps or overlaps of different regions can 
appear in the two-loop approximation. 

It is interesting to note that the scaling regimes that arise as solutions of the RG equations for the general model, 
in several cases correspond to interesting physical situations. In particular, the regimes governed by the points Q ,a 
correspond to the case of time-independent (or frozen) velocity field, while F ±,OL correspond to the rapid-change limit 
of the general model (2.5). 

To avoid possible misunderstandings we emphasize that the limits Uo — ► or go — ► 00, uq — > 00 are not supposed 
to be performed in the original correlation function (2.5); the parameters go, uo (and hence g, u) are fixed at some 
finite values. The behavior specific to the models (2.8), (2.9) arises asymptotically as a result of the solution of the 
RG equations, when the RG flow approaches the corresponding fixed point. This shows that in the regimes governed 
by the points Q^' a , the temporal fluctuations of the velocity field are asymptotically irrelevant in determining the 
inertial-range behavior of the passive field, which is then completely determined by the equal-time velocity statistics. 
In the regimes governed by R A , spatial and temporal fluctuations are both relevant, but the effective correlation time 
of the passive field becomes so large under renormalization that the correlation time of the velocity can be completely 
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neglected. The inertial-rangc behavior of the passive field is determined solely by the u> — mode of the velocity field; 
this is the case of the rapid-change model. In particular, this means that the coordinates of the fixed points and the 
anomalous exponents in such regimes must depend on the only exponent £ = 2e — rj or e that survives in the limit 
in question, and coincide with the corresponding dimensions obtained directly for the models (2.8) or (2.9). This is 
indeed the case, as one can see from the explicit expressions given above and in the next Section. 

As regards the value of the amplitude factor A in Eq. (2.1), it can remain an arbitrary parameter (for R A ), or 
can be attracted by one of the fixed points A* = ±1, or a (for Q or F); see Eq. (4.5). Again, these possibilities 
correspond to physical situations interesting as such. The case A* = 1 corresponds to the behavior characteristic 
of the magnetic model, where the pressure vanishes due to the relation (2.4) and the equation (2.1) becomes local. 
Therefore, the infrared behavior of the nonlocal model can be described by the same fixed point (or universality class) 
as that of the local (magnetic) model; the nonlocal pressure term does not affect the asymptotic properties of the 
passive field. The general model indeed becomes a "turbulence without pressure." The case A* = — 1 corresponds to 
the linearized NS equation with a given statistics of the background field (we recall, however, that this value of A* 
can be affected by the higher-order corrections in e and rf). Finally, the case A* — 0, a = corresponds the model, 
which (in its rapid-change variant) was introduced independently in a number of studies as an example of a linear 
system with pressure [53,37], a model with nontrivial mixing of composite operators [36,37] or a model which, with a 
proper choice of the forcing, can reproduce the anomalous exponents of the NS velocity field [54]. 

It is worth noting that for a ^ 0, the model with Ao = is not renormalizable, as follows from the analysis given 

(2) 

in Sec. III. That is, the second nonlinear term V { will be generated by the rcnormalization procedure, even if it was 
absent in the original equation (2.1). Of course, this fact does not mean that the model with Ao = and a ^ is 
inconsistent; it rather means that its IR behavior is described by one of the fixed point with A* ^ 0. 

V. CRITICAL SCALING. CRITICAL DIMENSIONS OF COMPOSITE OPERATORS 

Consider for definiteness some equal-time two-point quantity F(r) that depends on a single distance parameter 
r, for example, the pair correlation function of the primary fields 9, 9' or some composite operators. We assume 
that F(r) is multiplicatively renormalizable, i.e., F = ZpF R with certain renormalization constant Zp. Then the 
function F R (r) satisfies the RG equation of the form [Drq + jp] F( r ) = with the operator T>rq from (3.7) and 
jp = Vf+lnZp, cf. (3.8). The functions F and F R are equally suitable for studying the asymptotic behavior: the 
difference is in the normalization, choice of parameters (bare or renormalized) and the form of the perturbation theory 
(in g or m <?)■ The solution of the RG equation can be written in terms of the invariant variables introduced in (4.1). 
The analysis shows that in the IR asymptotic region, defined by the inequality Ar 3> 1 with A from (2.7) and any 
fixed mr with m = 1/L from (2.3), the invariant charges approach one of the IR attractive fixed points (the choice 
of the appropriate fixed point is discussed in the previous Section), and, as a result, the function F(r) takes on the 
self-similar form 

F(r) ~ nf A dF (Ar)- AF £(mr), (5.1) 

where d F and dp are the frequency and total canonical dimensions of F, respectively (see Sec. Ill and Table I), and 
£ is some function whose explicit form is not determined by the RG equation itself. The critical dimension Ap of the 
quantity F is given by the expression 

A F = dp + A aj dp+ip = d F - 7 *d£+ 7 *, (5.2) 

where -y F denotes the value of the anomalous dimension -fp at the fixed point in question, and A w = 2 — 7* with j K 
from (3.7) is the critical dimension of frequency. 

Each nontrivial fixed point Q ±,a , F ±,a and R A from Sec. IV corresponds to the scaling representation of the 
form (5.1) with its own set of critical dimensions Ap for all quantities F and A w . In general, these dimensions are 
infinite series in e and rj. For the rapid-change regime (that is, for the point R A ) they depend on the only exponent 
( = 2s — i] from (2.8), while for the regimes with quenched disorder (that is, for the points Q ±,a ) they depend on the 
only exponent e that survives in the limit (2.9). For the fixed point R A , the equation (3 X = determines the values 
of 7* = £ and A w = 2 — ( exactly, that is, without corrections of order ( 2 and higher. This follows from the explicit 
expressions (3.9), (4.12) and the vanishing of the anomalous dimensions 71.2; see the discussion below Eq. (4.12). For 
the fixed points F ± ' CI with w* ^ the equation (3 U — leads to the exact result 7* = 77; see Eq. (3.9). For the fixed 
point Q a with A* = a = 0, the exact result 7* = e follows from the equation fi g — and vanishing of 71; see the 
discussion in Sec. III. For the other regimes, the first-order expressions for 7* are directly obtained by substituting 
the coordinates of the fixed points (4.6a), (4.7a), (4.8a) into the one-loop expression (3.10c) for 7 K . The results can 
be summarized as follows: 
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7 K 



' ( = 2e — r] (exact) 
rj (exact) 
(d - 1 - a) 

£ (d-l) 
(d 2 -d- 3da + 2a + 2) 

8 (d-\)(d-2a) 

(d 2 - d + 2) 
£ d(d-l) 

(d 2 - da 3 + da 2 - a 3 + a 2 - a - 3) 

£ (d 2 - da 3 + da 2 + da - a 3 + a 2 + a - 3) 
k e (exact) 



for R A , 
for F ± - a , 

for Q+, 
for Q~, 

for Q~ and a = 0, 
for g Q , 

for Q Q and a = 0. 



(5.3) 



Then the dimensions Aj?, e.g., of the primary fields 9, 9' are obtained from Eq. (5.2) and the data from Table I, 

Aj = -l + 7 :/2, A 9 ,=d- 7 ;/2, (5.4) 

and the dimensions of their correlation functions are given by simple sums over the fields entering into the function. 

In the following, the crucial role will be played by the critical dimensions A[n, I] associated with the irreducible 
tensor composite fields ("local composite operators" in the field theoretic terminology) built solely of the fields 9 at 
a single spacetime point x — {t, x}. They have the form 



F[nJ]=6 n (x)---6 H (x) {6 i {x)e i {x)f + . . 



(5.5) 



where I < n is the number of the free vector indices and n = I + 2p is the total number of the fields 9 entering 
into the operator; the vector indices and the argument x of the symbol F[n, I] are omitted. The dots . . . stand for 
the appropriate subtractions involving the Kronecker delta symbols, which ensure that the resulting expressions are 
traceless with respect to contraction of any given pair of indices, for example, 9i9j —Sij9k0k/d, 9i9j9k — (5ij9k + Sik9j + 
5jk9i)9 2 1 (d + 2) and so on. We also note that the numbers n and I are even or odd simultaneously. 

Owing to the coincidence of the arguments, additional UV divergences arise in the correlation functions involving 
such operators. They are eliminated by means of the additional renormalization procedure, which gives rise to new 
(independent of ^1,2, « from (3.1)) renormalization constants. The analysis similar to that given in Refs. [41,43] for 
the scalar and in [33] for the magnetic model shows that, in the case at hand, these constants can be calculated in the 
model without forcing (the bare propagator (9i6j)o from Eq. (2.11) does not enter into the corresponding Fcynman 
diagrams). Then the operators (5.5) appear multiplicatively renormalizable: F[n, I] — Z[n, I] F R [n, I]. 

We have calculated all the constants Z[n,l] in the one-loop approximation (first order in g) in the MS scheme 
for the special case rj = and arbitrary e, which is sufficient to find the corresponding anomalous dimensions 
7 [n, /] = V^\nZ[n,l}; cf. the discussion above Eq. (3.4). We omit the calculation (which is very similar to that 
performed in Refs. [41,43,33] for the scalar and magnetic cases) and give only the final result: 



Z[n,l] = 1 + 



gS d 



8ed{d+2) (u+1) 



[.4 2 Qi + aQ 2 ] 



(5.6) 



where 



Qi = n(d + n)(d - 1) - 1(1 + d - 2)(d + 1), 
Q 2 = n(dn + n - d)(d - 1) - 1(1 + d - 2). 

For the anomalous dimension we thus obtain: 

According to Eq. (5.2), the critical dimensions of the operators (5.5) are given by 

A[n,l] = nA e + 7 *[n,Z] 



(5.7) 



(5.8) 



(5.9) 



with Ag from (5.4) and 7 *[n, I] is the value of the anomalous dimension (5.8) at one of the nontrivial fixed points 
(4.6a), (4.7a), (4.8a), (4.9a), (4.10a), (4.11a), (4.14a). Note that owing to the renormalization, the critical dimensions 
of the operators (5.5) differ from the naive sum of the dimensions of the fields 9 that constitute the operator. The 



17 



exception is provided by the case A* = a = 0, when A[n, Z] = nAg exactly (the proof is similar to that given, e.g., in 
Ref. [18] for the scalar case). 

The results for the anomalous dimensions 7* [n, I] can be summarized as follows: 

7>,Z] = -(2e-ri) [AlQ 1 +aQ 2 ] /2x 



x < 



(d + 2)(d-l + a) 
1 

{d - 2){d + 1 - a) 



(d 2 - da 3 + da 2 + 2da - a 3 + a 2 + 3a - 3) 
1 



I (d 2 - 3) 



for the regimes with finite correlation time, 



for F+, 
for F~, 
for F a , 

for F a and a = 



(5.10) 



7>,Z] = -£ [AlQ 1 + aQ 2 ] /2> 



x < 



(d + 2)(d-l) 
1 

(d- l)(d- 2a) 



(d 2 — da 3 + da 2 + da — a 3 + a 2 + a — 3) 
1 



I (d 2 - 3) 

for the regimes with quenched disorder and 

j*[n,l] = -(2e-r)) [A 2 Q 1 + aQ 2 ] /2x 



for Q+, 
for Q~, 
for Q a , 

for Q a and a = 



(5.11) 



x < 



1 



(d 2 + A 2 - aA 2 - adA 2 + adA + dA + da + 3a - 3) 
1 



I (d 2 + A 2 + dA - 3) 



for R A , 

for R A and a = 



(5.12) 



for the regimes with zero correlation time. We recall that A* takes on the values 1, —1 and a for the fixed points 
F and Q labelled by the superscripts +, — and a, respectively, while for the rapid-change regime A* — A remains 
an arbitrary parameter. We also note that the dimensions (5.11) depend on the only exponent e that survives in the 
limit (2.9), while the dimensions (5.12) depend on the only exponent £ = 2e — r\ that survives in the limit (2.8). The 
exponents (5.10) also depend only on £, which seems to be an artifact of the first-order approximation. We also note 
that the exponents (5.12) were derived earlier directly for the rapid-change model (2.8) for some special cases: A = 1 
and a = in [33] , A = 1 and arbitrary a > in [34] and arbitrary A and a — in [35] . 
From Eqs. (5.3), (5.4) and (5.9)-(5.12) it follows that 



A[n,l] = -n + 0(e). 



(5.13) 



Thus for all nontrivial fixed points, at least for small values of e ~ 77, one has A[n, I] < A[k, j] if n > k regardless of 
the relation between I and j. For a fixed value of n one has 



A[n, I] > A[nj] and j[n, I] > j[nj] if I > j, 



(5.14) 



as one can easily see from the original expression (5.8) for the anomalous dimension 7(71, 1], properties of the poly- 
nomials (5.7) and the fact that the combination g/(u + 1) = x/(w + 1) that enters into Eq. (5.8) is positive definite 
for all nontrivial fixed points. [We recall that the parameters x, w were introduced above Eq. (4.12) for the proper 
description of the rapid-change regimes.] Thus for fixed n, the dimension A[n, I] decreases monotonically with I and 
reaches its minimum for the minimal possible value of I, that is, I = if n is even and I = 1 if n is odd. 



18 



The hierarchy relations (5.13), (5.14) will be important in the discussion of the inertial-range behavior of various 
correlation functions, in particular, in the issue of the large-scale anisotropy persistence; see Sec. VI. Similar inequal- 
ities were established earlier for the case of the passive scalar field, advected by the velocity ensemble (2.5), (2.6) in 
Refs. [41,43], and for the magnetic field advected by the Kraichnan ensemble (2.8) in [39,33]. 

It also follows from (5.13) that the critical dimensions A[n, /] are negative and that the spectrum of their dimensions 
is not bounded from below; these properties, typical of the models of turbulence, will also be important in the following. 

If the random force / is introduced in Eq. (2.1), none of the formulas (5.6)-(5.12) change, because the bare correlator 
(99)o from Eq. (2.11), which becomes nonzero, does not enter the relevant diagrams. However, the operators (5.5) 
are no longer renormalized multiplicatively: the operator F[k,j] can admix in renormalization to the operator F[n, I] 
if, and only if, k < I. In general, j ^ I, but if the correlator (2.3) is isotropic, only operators with j = I can mix in 
renormalization. The admixture of operators with k > I is impossible due to the absence of appropriate diagrams; 
this is a consequence of the linearity of the original equation (2.1) in 9 and /. The admixture of operators with k = n 
and j ^ I is also impossible, because the corresponding diagrams do not involve the correlator (69) q and therefore do 
not "feel" the violation of the rotational symmetry caused by the function (2.3). 

As a result of the mixing, the operator F[n, I] becomes a finite sum of contributions with definite critical dimensions 
A[n, 1} and A[fc,j] with k < n and, in general, all possible values of j allowed for a given k. However, due to the 
relation (5.13), the leading term is still given by the original contribution with the dimension A[n, I], while the new 
contributions with k < n give only corrections that vanish in the IR range Ar 3> 1 in expressions like (5.1). In what 
follows, we shall be interested only in the leading terms of and thus we can ignore the mixing and treat the operator 
F[n, I] as if it has the definite critical dimension A[n, I]. 

VI. OPERATOR PRODUCT EXPANSION AND THE ANOMALOUS SCALING 

The representation (5.1) for any scaling function £(mr) describes the behavior of the correlation function F(r) for 
Ar 3> 1 and any fixed value of mr. The inertial range corresponds to the additional condition that mr <C 1. The 
form of the function ((mr) is not determined by the RG equations themselves; in the theory of critical phenomena, 
its behavior for mr — > is studied using the well-known Wilson operator product expansion (OPE); see, e.g., Ref. [7]. 
This technique is also applicable in the theory of turbulence; see, e.g., Ref. [23-25]. 

According to the OPE, the equal-time product Fi(x)F 2 (x') of two renormalized composite operators at x = (x + 
x')/2 = const and r = x — x' — > can be represented in the form 

F 1 (x)F 2 (x') = ^C F (r)F(i,x), (6.1) 

F 

where the functions Cp are the Wilson coefficients regular in m 2 and F are, in general, all possible renormalized 
local composite operators allowed by symmetry; more precisely, the operators entering into the OPE are those which 
appear in the corresponding Taylor expansions, and also all possible operators that admix to them in renormalization. 
If these operators have additional vector indices, they are contracted with the corresponding indices of the coefficients 

Cp. 

Without loss of generality it can be assumed that the expansion in Eq. (6.1) is made in the operators with definite 
critical dimensions A^. The renormalized correlation function (F\(x)F2(x')) is obtained by averaging Eq. (6.1) with 
the weight expS^ with Sr from Eq. (3.1), the quantities (F) appear on the right hand side. Their asymptotic 
behavior for m — > is found from the corresponding RG equations and has the form (F) oc m Ap . 

From the operator product expansion (6.1) we therefore find the following expression for the scaling function ((mr) 
in the representation (5.1) for the correlation function (Fi(x)F2(x')): 

Z(mr) = Y, A F( mr ) AF ' ( 6 - 2 ) 

F 

where the coefficients Ap(mr) arc regular in (mr) 2 . 

The quantities of interest are, in particular, the equal-time pair correlation functions of the composite operators 
(5.5). For these, the representation (5.1) is valid with the dimensions d F = — (n + k)/2, dp = — (n + k) and 
Ap = A[n,l] + A[k,j] with A[n,l] from (5.9): 

(F[n,i](t,x)F[ft,j](t,x')) = (K )- {n+k)/2 K- {n+k) (kr)-^ n ^- A ^Z{mr) (6.3) 
where Ar » 1 and £ (mr) is the corresponding scaling function. 
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As already said above, the operators entering into the OPE are those which appear in the corresponding Taylor 
expansions, and also all possible operators that admix to them in renormalization. The leading term of the Taylor 
expansion for the function (6.3) is given by the n-th rank tensor F[n + k, I + j] from Eq. (5.5). Its decomposition in 
irreducible tensors gives rise to the operators F[n + k,p] with all possible values of p < I + j; the admixture of junior 
operators (see the end of Sec. V) gives rise to all the monomials F[s,p] with s < n + k and all possible p allowed for 
a given s. Hence, the asymptotic expression for the structure function £(mr) for mr <C 1 has the form 

n+k s 

£(mr) = E K> (™0 A[S ' P] +■■■], 

s— p—Ps 

with the dimensions A[fc,p] from Eq. (5.9). Here and below p s denotes the minimal possible value of p for given s, 
i.e., p s = for k even and p s = 1 for k odd; A sp are some numerical coefficients dependent on the parameters like 
e, d and so on. The dots in Eq. (6.4) stand for the contributions which arise from the composite operators that, in 
addition to the field 6, involve the other fields 8' , v and/or derivatives d, dt- 

The leading term of the expression (6.4) for mr 1 is determined, obviously, by the minimal possible dimension 
Ap that appear on its right-hand side, provided this minimal dimension exists. In our model, there are infinitely 
many operators with negative critical dimensions, and the spectrum of their dimensions is not bounded from below. 
(It is possible to show on general grounds that, if a model involves one negative dimension, it necessarily involves 
infinitely many negative dimensions with unbounded spectrum.) If all these operators appeared on the right-hand 
side of the representation (6.4), we would have to sum up their contributions in order to find the asymptotic behavior 
at mr — > 0. This problem is indeed encountered for the stochastic NS equation [23], and is discussed in Refs. [25,24] 
in detail. 

In our model, however, there is no such problem, at least for small e. The contributions of the operators F[s,p] 
with s > n + k (which would be more important) do not appear in Eq. (6.4), because they are absent in the Taylor 
expansion of the correlator (6.3) and do not admix in renormalization to the terms of the Taylor expansion; see 
Sec. V. As already noted there, this is a manifestation of the linearity of the original equation (2.1) in 8 and /. 
What is more, one can show that for any operator F that appear in the OPE (and not only the operators (5.5) built 
solely of the fields 8) the number of the fields 8 cannot exceed the total number of the fields 8 on the left-hand side; 
therefore their dimensions cannot appear in (6.4). It then follows that the leading term in (6.4) is determined by an 
operator built solely of the fields 8 and containing the maximal possible number of the fields, that is, n + k. The 
operators containing less than n + k fields 8 give only corrections, as follows from the hierarchy relation (5.13). The 
operators involving the fields 8' , v and/or derivatives also give only corrections, because the canonical dimensions 
of these additional factors are positive (see Table I) and thus increase the total canonical dimension of the operator 
in comparison with the corresponding operator built solely of the fields 8. Furthermore, from the hierarchy relation 
(5.14) it follows that, for the fixed number of the fields 8, the minimum of the dimension is achieved for the minimal 
number of vector indices, that is, for the scalar operator F[n + k, 0] if the sum n + k is even and the vector operator 
F[n + k, 1] if the sum n + k is odd. 

We therefore conclude that the leading term of the small-mr behavior of the scaling function (6.4) has the form 
£ ~ [mr)^ n+k ' ln+k ^. Substituting this expression into Eq. (6.3) gives the desired leading term of the correlation 
function of two operators (5.5) in the inertial-range (Ar 3> 1, mr <C 1): 

(F[n,l](t,y:)F[k,j}{t,x')) ~ ( Ko )-(™+W2 A -(«+fc) (Ar) -A[rM]-A[M ^Aln+k,l n+k ] (g 5) 

with the dimensions A[n, I] from (5.9). 

VII. ANOMALOUS SCALING IN ANISOTROPIC SECTORS. HIERARCHY OF ANISOTROPIC 

CONTRIBUTIONS 

Without loss of generality, it can always be assumed that the expansion (6.1) is made in irreducible traceless tensor 
composite operators. Then averaging Eq. (6.1) with the weight expSu automatically produces the decomposition of 
the correlation function in irreducible representations of the rotation group SO(d), similar to that employed e.g. in 
Refs. [57-62] for description of the NS turbulence. This becomes especially clear if the left-hand side of Eq. (6.1) 
involves only scalar quantities and the anisotropy, introduced by the correlator (2.3), is uniaxial, that is, specified 
by a single constant unit vector n. Then the mean value (F) of a Z-th rank tensor operator F on the right-hand 
side of (6.1) is an irreducible traceless Z-th rank tensor built only of the vector n and the Kronecker delta symbols. 
Its vector indices are contracted with the indices of the corresponding Wilson coefficient Cf(t), which gives rise to 
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(6.4) 



the Legendre (or Gegenbauer for arbitrary d) polynomial of order I. In general, decomposition in (hyper)spherical 
harmonics (see e.g. [64] and references therein) or its analogs for tensor quantities (see e.g. [65] and the references) 
will be encountered. 

The rank I of the operator can be viewed as the measure of anisotropy of the corresponding contribution in expansion 
(6.4). If the forcing is isotropic, that is, the function C(r) in the correlator (2.3) depends only on r = |r|, only scalar 
operators with / = have nonvanishing mean values, and only their dimensions appear on the right-hand side of Eq. 
(6.4). In general, tensor operators with / ^ also contribute to (6.4). Owing to the relations (5.14), the leading term 
of the asymptotic behavior at mr — > is still given by the scalar operator with Z = (it has the minimal dimension 
among the operators with a fixed number of the fields). We thus conclude that the leading term is given by the same 
expression (6.5) for both the isotropic and anisotropic forcing, while anisotropic contributions with I > give only 
subleading terms (corrections). What is more, relations (5.14) show that these contributions reveal a kind of hierarchy 
related to the degree of anisotropy: the higher is the rank of the operator, the less important is its contribution to 
the inertial-range behavior. 

For the first time, the hierarchy relations for anisotropic contributions were derived in Ref. [39] for the magnetic 
field, passively advected by Kraichnan's velocity ensemble (2.8), and in Ref. [41] for the scalar field, advected by the 
Gaussian velocity field specified by the correlator (2.6) (in both cases with general d and a = 0). In the first of these 
papers, anomalous exponents were found exactly for the pair correlation function, while in the second the exponents 
were derived only in the one-loop approximation, but for all the higher-order correlation functions. Later these results 
were reproduced in Ref. [40] for the magnetic field (only for d = 3, but also including helical contributions) and in 
[52] for the scalar field advected by Kraichnan's ensemble. Generalization to the higher-order correlation functions of 
the magnetic field was given in [33] (for a = 0), while generalizations to the case of general a were obtained in [34] 
(magnetic field and Kraichnan's ensemble) and in [43] (scalar field and the ensemble (2.6)). Generalization to the 
general vector model (2.1) and Kraichnan's ensemble was given in [35]. 

So far, analytical results of such kind have been obtained only for passive fields, advected by the synthetic Gaussian 
velocity ensembles. However, numerical simulations and real experiments show that the picture outlined above 
appears rather general, being observed by the passive scalar field advected by the two-dimensional NS field in the 
inverse energy cascade [63] and by the NS velocity field itself [59-62]. These observations justify and make more precise 
old phenomenological ideas about the isotropization of the inertial-range turbulence in the presence of a large-scale 
anisotropy. Nevertheless, the anisotropy survives in the inertial range and reveals itself in odd correlation functions, 
in disagreement with what was expected on the basis of the cascade ideas. We shall return to this important issue 
in the Conclusion, and now let us briefly discuss the influence of compressibility on the hierarchy of the anomalous 
exponents. 

Effects of the compressibility on the anomalous scaling in anisotropic sectors were studied earlier for the scalar [31] 
and magnetic [34] fields advected by Kraichnan's ensemble (2.8) and for the scalar advected by the velocity ensemble 
(2.6) with finite correlation time [43]; see also Ref. [66] for a summary. In all those cases the conclusion was the same: 
the hierarchy expressed by the relation (5.14), which can be rewritten as 

dA[n,p]/dp > 0, (7.1) 

remains valid for all values of the compressibility parameter < a < +oo, but it always becomes less pronounced as 
a grows, 

d 2 A[n,p]/dpda < 0. (7.2) 

This means, in particular, that the anisotropic corrections in Eq. (6.4) become closer to each other and to the leading 
term as a grows. Thus the compressibility enhances the penetration of the large-scale anisotropy into the inertial 
range. This penetration is even more manifest for the odd-order ratios of the correlation functions: the skewness factor 
grows for mr — ► 0, provided a is large enough, while the growth of the hyperskewness factor and other higher-order 
ratios becomes much faster than for the incompressible case; see the discussion in [31,43,34,66]. 

No such definite conclusions can be drawn for the general vector model. The straightforward analysis of the explicit 
expressions (5.10)-(5.12) shows that the derivative d 2 A[n 7 p]/dpda is negative and, therefore, the behavior described 
above takes place only in two regimes described by the fixed points F + (always) and R A (only if the relation 

-3 + Ad + d 2 - A 3 d(d+1) + A 4 (d+1) 2 - A 2 (d 2 + Ad+2) > 0, (7.3) 

which is independent on a, is satisfied). For all the other cases, one finds d 2 A[n,p]/dpda > and the behavior is 
opposite: compressibility suppresses the penetration of the large-scale anisotropy into the inertial range, anisotropic 
contributions become further from one another and from the isotropic term. 



21 



It is tempting to attribute this "inverse behavior" to the combined influence of the compressibility and pressure. 
Indeed, the "normal" (scalar-like) behavior takes place for the magnetic regime F + , for which A* = Ao = 1 and the 
pressure term (2.4) vanishes. For the rapid-change regime i?" 4 and reasonable values of d, the relation (7.3) is satisfied 
only in the restricted area around the point A* = 1 (including A* = 1, in agreement with the analysis of [34]), where 
the pressure effects arc relatively small. However, the magnetic regime in a frozen velocity field, Q + , demonstrates 
the inverse behavior. 

The dependence on the parameter A, that controls the pressure effects, is essentially different from the dependence 
on a: for the regimes with nonzero (finite or infinite) correlation time, the value of the corresponding invariant variable 
can take only three discrete values A* = — 1, 1, a. The value of A* which is realized for a given regime depends 
on a but not on A; see Eq. (4.5). The case A* = 1 corresponds to the magnetic (pressureless) equations; then the 
general model indeed becomes a "turbulence without pressure," despite the presence of the nonlocal pressure term 
in the original stochastic equation. For the zero-correlated regime i?- 4 , the anomalous exponents retain a continuous 
dependence on A; for the incompressible case (a = 0) it was discussed in [35]. The value of A — 1, where the 
pressure effects disappear, is not distinguished at all; the derivative d 2 A[n,p]/dpdA at A = 1 is positive for almost 
all parameters (namely, for d > — 1.3a + 1.5, the approximate relation obtained numerically), but is negative definite 
e.g. for A = 0. 

VIII. CONCLUSION 

We have studied a model of a divergence-free (transverse) vector quantity passively advected by a random Gaussian 
velocity field with finite (and not small) correlation time. The model is described by an advection-diffusion equation 
with a random large-scale stirring force, nonlocal pressure term and the most general form of the incrtial nonlinearity. 
The correlation function of the advecting field mimics some properties of the real inertial-range turbulence: the energy 
spectrum has the form E{k) oc k x ~ 2e , while the correlation time scales as k~ 2+v . An advantage of the model is the 
possibility to control the pressure contribution and thus study its effects on the inertial-range behavior. Another 
reason to study the general case is the possibility to describe in a uniform way several special cases interesting on 
their own: the kinematic magnetic model, linearized NS equation and the special model without the stretching term, 
which possesses additional symmetry and has a close formal resemblance with the nonlinear NS case. 

We have shown that the system exhibits various types of inertial-range asymptotic behavior, characterized by 
nontrivial anomalous exponents; the latter are analytically calculated to first order in e ~ rj, including the anisotropic 
sectors. 

The key points of our analysis are the existence of a field theoretic formulation of the original stochastic problem 
(Sec. II), multiplicative renormalizability of the corresponding field theory (Sec. Ill), existence of nontrivial IR- 
attractive fixed points of the corresponding RG equations in the physical region of the parameters (Sec. IV) and the 
possibility to identify the anomalous exponents with the critical dimensions of certain composite operators (Sees. V 
and VI) . This allows one to construct a systematic perturbation expansion for the exponents; the practical calculations 
have been performed to the first nontrivial order in e ~ rj (one- loop approximation). 

Existence of explicit one-loop expression allows one to discuss the stability scaling regimes and the universality of 
the corresponding exponents, that is, their (in)dependence on the pressure, anisotropy, compressibility, forcing and 
so on, or, more technically, on the exponents e, rj and the amplitudes Ao, go, uq and a in the stochastic equation 
(2.1) and the correlator (2.6) of the advecting velocity. Although the behavior of the vector model is much richer 
than that of its scalar counterpart, the general picture appears essentially the same: the exponents are universal in 
the sense that they depend on the exponents e and rj, but do not depend on the amplitudes in (2.6) and the forcing 
(2.3); the exponents related to anisotropic contributions show a hierarchy related to the degree of anisotropy (more 
anisotropic contributions are less important); this hierarchy holds for all scaling regimes, regardless of the values of 
the compressibility parameter a from (2.5) and the pressure parameter A from (2.1). Consider these points in more 
detail. 

Scaling regimes and universality classes. Infrared asymptotic behavior of our model is completely described 
by seven different scaling regimes, or universality classes, each corresponding to a set of anomalous exponents. For 
the given set of the parameters e, r\ and a, only one of these regimes can be realized, irrespective of the values of 
the amplitudes Ao, go, uo- Three regimes correspond to finite correlation time of the advecting velocity field, and 
three regimes correspond to infinite correlation time (or time- independent velocity). Each of these two sets involve 
the magnetic (pressureless) case, linearized NS equation and the model with a-dependent effective amplitude in front 
of the stretching term; for a = the latter gives the special model with no stretching term, whose rapid-change 
version was studied e.g. in [37,53,54]. The scaling exponents for these regimes depend on the exponents e, rj and the 
amplitude a, but they are independent of the values of the amplitudes Ao, go, uq. The remaining seventh regime 
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corresponds to the rapid-change velocity (zero correlation time), the corresponding exponents depend also on Ao- 

To avoid possible confusion we stress that the behavior specific to the aforementioned classes, e.g. magnetic model 
with infinite correlation time, arises automatically when the RG flow approaches the fixed point which is IR attractive 
for the given choice of parameters e, r\, a; the frozen limit (2.9) or the substitution Ao = 1 arc not performed 
in the original model and the parameters Ao, go, u are fixed at such finite values. In particular this means that 
the anomalous exponents in those regimes are independent of the correlation time (more precisely, the ratio of the 
correlation time of the velocity field and the turnover time for the scalar field, measured by the parameter uo; see the 
discussion in [41]). In this sense, one can speak about the universality of the anomalous exponents in our model. 

Independence of the forcing. Zero-mode picture. As we have seen, the critical dimensions of all composite 
operators (5.9), and therefore the corresponding anomalous exponents (including anisotropic sectors), are independent 
of the forcing, specified by the correlator (2.3). In particular, this means that they remain unchanged, when the stirring 
force in Eq. (2.1) is replaced by the imposed mean constant field, like in Refs. [33,39]. The role of the forcing is to 
maintain the steady state of the system and thus to provide nonzero amplitudes for the power-like terms with those 
universal exponents. 

This behavior is already well known for the passive scalar fields [41,43] advected by the velocity (2.6) or vector 
fields, advected by the zero-correlated velocity [33,39]. 

In the language of the RG (which is equally applicable to the case of a zero or finite correlation time) this is explained 
as follows: the stirring force or the mean field do not enter into the diagrams that determine the renormalization of the 
operators (5.9), so that their dimensions are independent of the forcing. Similar diagrams determine the contributions 
of those operators into the operator-product expansions (6.1), which are nontrivial even for the unforced model. The 
difference is that for the unforced model, mean values of the operators vanish, and they give no contribution to 
the right-hand sides of representations like (6.4). For the isotropic correlator (2.3), scalar operators acquire nonzero 
mean values and contribute to the right-hand side of (6.4), while for the anisotropic correlator or the imposed mean 
field, the mean values of irreducible tensor operators also become nonzero and their contributions are "activated" in 
representations (6.4). 

For the case of zero correlation time, when the equal-time correlations functions satisfy exact closed differential 
equations, the above picture it is easily understood in the language of the zero- mode approach [14]: forcing terms 
do not affect the corresponding differential operators; thus the anomalous exponents, determined by the zero modes 
(solutions of homogeneous unforced equations) also are independent of the forcing. On the contrary, the amplitudes 
are determined by the matching of the inertial-range zero-mode solution with the forced large-scale solutions, which 
is only possible in the presence of the forcing terms. 

The exact resemblance in the behavior of the rapid-change models and the finite-correlated cases suggests that 
for the latter, the concept of zero modes (and thus of statistical conservation laws) is also applicable, although the 
corresponding equations are not differential and involve infinite diagrammatic series. 

Hierarchy of anisotropic contributions. In the presence of the large-scale anisotropy (that is, the anisotropy 
introduced at scales of order L by the forcing) , correlation functions of the model can be decomposed in irreducible 
representations of the ci-dimensional rotation group SO(d). Such a decomposition naturally arises from the corre- 
sponding OPE, provided it is made in irreducible traceless tensor composite operators; the rank I of a tensor operator 
can be used to label the terms of the SO (d)-expansion and can be viewed as the measure of anisotropy of the cor- 
responding term ("sector"). Thus each anisotropic sector is characterized by its own set of scaling exponents, the 
leading term is given by the l-th rank composite operator with minimal critical dimension. 

Explicit expressions for these dimensions were obtained to the first order in e and rj. They reveal an hierarchy 
related to the degree of anisotropy: the higher is the rank of the operator (the more anisotropic is the contribution), 
the larger is the corresponding dimension, and thus the less important is its contribution to the inertial-range behavior. 

This hierarchy, expressed by the relations (5.14) or (7.1), holds for all nontrivial scaling regimes of our model, all 
values of the parameters a, A, d and so on. It is similar to the hierarchy relations derived earlier for the passive scalar 
[41,43,52] and magnetic fields [33,39,40] advected by the Gaussian velocity ensembles. 

In particular, this means that the overall leading term is given by the exponent from the isotropic sector, ant it 
is therefore the same for the isotropic and anisotropic forcing. It also should be stressed that the independence of 
the scaling behavior in different sectors is a direct consequence of the linearity of our model, independence of the 
exponents on the random force, and the SO{d) symmetry of the unforced model. On the contrary, the hierarchy of 
the exponents follows from the explicit expressions, obtained only by practical calculation. 

According to the Kolmogorov-Obukhov theory [1,2], the anisotropy introduced at large scales by the forcing (bound- 
ary conditions, geometry of an obstacle etc) dies out when the energy is transferred down to smaller scales owing to the 
cascade mechanism (isotropization of the developed turbulence in the inertial-range). The analytical results discussed 
above confirm this classical concept and give a more quantitative picture of the isotropization. The relevance of these 
results for more realistic situations (scalar advected by the two-dimensional NS field or the turbulent velocity itself) 
is briefly discussed below. 
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Effects of compressibility. The anomalous exponents explicitly depend on the parameter a > that measures 
the compressibility of the fluid. For the regimes determined by the fixed points F + (magnetic model with finite 
correlation time) and i?" 4 (zero correlation time, with additional inequalities for the parameter A satisfied by the 
magnetic case A = 1 and its vicinity), the hierarchy of anisotropic contributions becomes less pronounced as a grows: 
the anisotropic corrections in Eq. (6.4) become closer to each other and to the leading term as a grows. Thus the 
compressibility enhances the penetration of the large-scale anisotropy into the inertial range. The situation is opposite 
for all the other regimes, which arguably can be attributed to the influence of the pressure term. 

Effects of pressure. The dependence on the parameter A, that controls the pressure effects, is essentially different 
from the dependence on a: for the regimes with nonzero correlation time, the value of the corresponding invariant 
variable can take only discrete values A* = —1, 1, a. The behavior for the fixed point A* = 1, which corresponds to 
the magnetic case ("turbulence without pressure"), shows no serious difference from the regimes with pressure. For 
the rapid-change limit, the exponents continuously depend on A, and the value of A = 1, where the pressure effects 
vanish, is not distinguished either. 

Relevance for the NS turbulence. The picture outlined above for passively advected fields (a superposition 
of power laws with universal exponents and nonuniversal amplitudes) seems rather general, being compatible with 
that established recently in the field of NS turbulence, on the basis of numerical simulations of channel flows and 
experiments in the atmospheric surface layer; see Rcfs. [57-62] and references therein. It was shown that the leading 
terms of the inertial-range behavior are the same for isotropic and anisotropic forcing [57,58]. In the papers [59-62], 
the velocity correlation functions were decomposed in the irreducible representations of the rotation group. It was 
argued that in each sector of the decomposition, scaling behavior can be found with apparently universal exponents. 
The amplitudes of the various contributions are nonuniversal, through the dependence on the position in the flow, 
the local degree of anisotropy and inhomogeneity, and so on. 

This is rather surprising because the equations for the correlation functions in such cases are neither closed nor 
isotropic and homogeneous. Although the hierarchy similar to Eq. (5.14) is demonstrated by the critical dimensions 
of certain tensor operators in the stirred NS turbulence, see Sec. 2.3 of [25], the relationship between them and 
the anomalous exponents is not obvious there. It is worth recalling here that the so-called "additive fusion rules," 
hypothesized for the NS turbulence in a number of papers, Refs. [16,17,26], and characteristic of the models with 
multifractal behavior (see Ref. [27]), arise naturally in the context of the models of passive advection owing to their 
linearity. The existing results for the Burgers turbulence can also be interpreted naturally as a consequence of similar 
fusion rules, where only finite number of dangerous operators contributes to each structure function, see Ref. [28]. 

One can thus speculate that the anomalous scaling for the genuine turbulence can also appear a linear phenomenon 
in the following sense. Let us split the total velocity field into the two parts, the background field and the perturbation 
(e.g., large-scale and small-scale, or soft and hard components), linearize the original stochastic equation with respect 
to the latter, choose an appropriate statistics for the former (e.g. Gaussian distribution with Kolmogorov exponents, 
the description suggested for the large-scale field by the experiment). Then the small-scale perturbation field will 
show anomalous scaling behavior with nontrivial exponents, which can be calculated systematically within a kind of 
e expansion. The corrections due to the nonlinearity can be treated perturbatively, and if they appear irrelevant (e.g. 
in the sense of Wilson), they will not affect the exponents calculated within the linearized model. In such a case the 
passive vector field can give the anomalous exponents for the NS velocity field exactly. In other words, such linearized 
model will belong to the same universality class as the real NS equation, like the simplified Ising or Heisenberg models 
are believed to belong to the same universality class as real ferromagnets or binary alloys. It thus might happen that 
the anomalous behavior of the real inertial-range turbulence is exactly described by one of the nontrivial fixed points 
for the passive vector model. 

Of course, one should not insist too much on such a simple scenario for the anomalous scaling, but it is worthy 
of attention. In this connection, we could also recall that the passive vector field can indeed reveal the anomalous 
exponents of the stochastic NS velocity field if the random forcing of the former is chosen to be statistically correlated 
with that of the latter; see [54]. 

Validity of the e expansion and the applicability of the model. A serious question is that of the validity 
of the e expansion and the possibility of the extrapolation of the results, obtained within the e expansions, to the 
finite values e = 0(1). For the rapid-change model, the e expansion works surprisingly well. It was shown [29] that 
the knowledge of three terms allows one to obtain reasonable predictions for finite e ~ 1; even the plain e expansion 
captures some subtle qualitative features of the anomalous exponents established in analytical and numerical solutions 
of the exact zero-mode equations and numerical experiments. The agreement can be further improved by using special 
tricks (like the "inverse" e expansion) or interpolation formulas [29]. 

In the case of the Gaussian model with a finite correlation time, however, there is a natural upper bound for the 
range of validity of the results, obtained within the e expansion: for e > 1 the velocity field (and hence all its powers) 
become dangerous (its critical dimension A v — 1 — <s, known exactly due to the Gaussianity, becomes negative). The 
spectrum of their dimensions is unbounded from below, and in order to find the small-mr behavior one has to sum up 
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all their contributions in the representations like (6.4). This problem is discussed in detail in [41] for the passive scalar 
field; the infrared perturbation theory was employed there to perform the required summation for the pair correlation 
function, in the frozen regime, and within the one-loop approximation for the Wilson coefficients. It was argued that, 
in that special case, anomalous behavior is described by the same exponent below and above the boundary e = 1, but 
in general the problem remains open. 

Physically, this is a manifestation of the fact that for e > 1, the so-called sweeping effects (kinematic transfer of 
the small-scale turbulent eddies by the large-scale ones) become important. In a Galilean-covariant problem such 
composite operators would not give any contribution into the Galilean invariant quantities (structure functions), as it 
happens in the RG approach to the stochastic NS equation; see the discussion in Refs. [23] and [45]. As was pointed 
out in Rcf. [46], the Gaussian model with finite correlation time suffers from the lack of Galilean invariance and 
therefore misrepresents the sweeping effects: they penetrate into the correlation functions of the scalar and can lead 
to their strong unphysical dependence on L. Therefore the value s = 1 can also be viewed as the threshold above 
which the model itself becomes unphysical. [To justify the Gaussian model for e > 1, however, one may recall that 
the results of [46] show that it gives a reasonable description of the passive advection in an appropriate frame, where 
the mean velocity field vanishes.] 

We may therefore conclude that the next important step is the analytical derivation of anomalous exponents of a 
passive scalar and vector quantities advected by the Galilean covariant velocity ensemble, generated by the stochastic 
NS equation; this work is now in progress. 
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FIG. 1. Regions of infrared stability of the three fixed points (4.6a), (4.7a), and (4.8a) corresponding to quenched disorder 
for a representative set of values of the space dimension d and the parameter a of the relative strength of the longitudinal part 
in the correlation function of the velocity field. The region of stability for any indicated values of d and a lies between the 
dashed ray (e = 0, ij < 0) and the correspondingly marked dash-dotted ray in the upper half of the 77, e plane. Note different 
scales on the coordinate axes. 
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FIG. 2. Regions of infrared stability of the fixed point (4.9a) with a finite correlation time in the asymptotic regime. Basins of 
attraction are shown for a representative set of values of the space dimension d and the parameter a of the relative contribution 
of the longitudinal part in the correlation function of the velocity field. The region of stability for any indicated values of d 
and a lies between the dashed ray (e = 77, ij > 0) and the correspondingly marked dash-dotted ray in the upper half of the r/, 
e plane. Note different scales on the coordinate axes. 
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FIG. 3. Regions of infrared stability of the fixed point (4.10a) with a finite correlation time in the asymptotic regime. Basins 
of attraction of the three fixed points are shown for a representative set of the longitudinal parameter a in space dimensions 
two, three and four. The region of stability for any indicated value of a lies from the dashed ray (e = r), r\ > 0) to the left up to 
the correspondingly marked dash-dotted ray. Note different scales on the coordinate axes. Contrary to the quenched-disorder 
case, negative values of the correlation falloff parameter e are (formally) allowed. 
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FIG. 4. Regions of infrared stability of the fixed point (4.11a) with a finite correlation time in the asymptotic regime. Basins 
of attraction of the three fixed points are shown for a representative set of the longitudinal parameter a in space dimensions 
two, three and four. The region of stability for any indicated value of a lies from the dashed ray (e = r), rj > 0) to the left up 
to the correspondingly marked dash-dotted ray. Note different scales on the coordinate axes. Negative values of e appear here 
as well as for the fixed point (4.11a). 



29 



